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MULTI-SOLITONS AND RELATED SOLUTIONS 
FOR THE WATER-WAVES SYSTEM 

MEI MING, FREDERIC ROUSSET, AND NIKOLAY TZVETKOV 

Abstract. The main result of this work is the construction of multi-solitons solutions that is to 
say solutions that are time asymptotics to a sum of decoupling solitary waves for the full water 
(.^ ■ waves system with surface tension. 

(N' 

1. Introduction 
00 

We consider the motion of an irrotational, incompressible fluid with constant density. We consider 
the situation where the fluid domain is a strip with a rigid bottom and a free surface: 

n, ; 

^ . Oj = {y = (X, z) G M'^+i : -H <z< r]{t, X)}, 

,^ I where t is the time, d = 1,2 is the horizontal dimension, if is a parameter defining the fixed 

C^ ' bottom z = —H and z = r]{t, X) is the equation of the unknown free surface at time t. We 

shall say that we are in the one dimensional case when X = x £ M and in the two-dimensional 
case when X = {x,y) £ M?. A large part of the paper will be devoted to the one-dimensional 
situation. We denote by u the speed of the fluid, since the motion is irrotational, it is given by 
^ ! u = Vy <1> = (Vx^, dz^) for some scalar function $ and hence we find that inside the fluid domain 

(N| : Vyn = Ay$ = (Ax + 9^)^ = 0. (1.1) 
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On the boundaries of 0,t, we make the usual assumption that no fluid particles cross the boundary. 
At the bottom of the fluid this reads 

dz'^{t,X,-H) = (1.2) 

and on the free surface, this yields the kinematic condition 

dt7]{t, X) + VxHt, X, 7]{t, X)) • VxTlit, X) - dMt, X, r,{t, X))=G. (1.3) 

On the free surface, we also need to impose the pressure, taking into account the surface tension 
and using the Bernouilli law to eliminate the pressure, we find that: 

dMt,X,r,{t,X)) + i|Vy$(t,X,,?(t,X))|2 +5r?(t,X) = hV x ■ . ^fj^^'f^^,, ^ • (1-4) 

The number b is the surface tension coefficient and g is the gravitational constant. The term 
gr]{t, X) is the trace of the gravitational force gz on the free surface. 

It is classical to rewrite the system (jl.ip . (jl.3p . ()1.4p as a system involving unknowns defined on 
the free surface only [31]. For that purpose, let us define the following Dirichlet-Neumann operator: 
for given r]{X) ^{X), we define $(X, z) as the (well-defined) solution of the elliptic boundary value 
problem 

(Ax + 92)c^ = 0, in {{X,z) : -H<z<ri{X)], <^{X,j^{X)) = ^{X), dMX,-H) = 0, 



and we define the Dirichlet-Neuniann operator as 

= Vl + |Vx??p(Vx,.^-n)U=^(x), 

where n is the unit outward normal vector on the free surface at the point z = ri{X). 
This allows to rewrite the system only in terms of the unknowns 

(r?(t, X), ^{t, X)) := {r,{t, X), ^t, X, r,{t, X))) . 

In the one-dimensional case, the ID water-wave problem can thus be written as 

dtrj = G[r]]ip 

By introducing the notations U = (rj, (/?)* and 



we shall write the water-wave system (jl.Sp in the abstract form 

dtU = T{U). (1.6) 

We know from [^ that for suitable parameters g, b and h, there exist solitary wave solutions 
Qc{x — ct) = {rjc{x — ct), (pc{x — ct)Y at speed c ~ y/gH. Here is a precise statement. 

Theorem 1.1 (Amick-Kirchgassner [3]). Suppose that 

„ = ^ = l+e^ ,3 = ^>^. (1.7) 

Then there exists Eq such that for every e £ (0, Eq) (which fixes the speed) there is a solution of 
il.5\) under the form 

Qc{x - ct) = {r]c{x - ct), ipc{x - ct)y = (^Hrje{H-\x - ct)), cHipe{H-\x - ct))) 

with 

r]e{x) = e^Qi{ex,e), ips{x) = £92(ex,e), 
where 0i and 02 satisfy: 

3d>0, Va>0, 3Ca>0, V (x,e) e M x (0,eo), |(9°ei)(x,e)| < C^e"'^!^! 
and 

3d>0, Va>l, 3Ca>0, V (x,e) G R x (0,eo), |(9^02)(2;,e)| < C^e""^!^! . 
Moreover Gi is even and G2 is odd. 

The main aim of this paper is to construct multi-solitons type solutions for the water-waves 
system (jl.Sp . This means that we want to construct a solution of p.Sp that tends to a sum of 
solitary waves with different speeds when t goes to infinity. For the sake of readability of the 
paper, we shall only consider the case of two solitary waves, the extension to an arbitrary numbers 
is straightforward (our proof will not use any particular symmetry or specificity related to the 
2-solitons case) . The construction of such solutions for semi-linear equations like the KdV equation 
or the Nonlinear- Schrodinger equation has been intensively studied recently, we refer for example 
to |18l [T9| [T6| [9l El [22] . We also refer to an earlier closely related work by Merle [20] which seems 
to initiate this line of research. 



In this paper, we shall show that the construction can be also performed for quasilinear equations 
(or even fully nonlinear equations) by focusing on the physically interesting example of the water- 
waves system (|1.5p . The new difficulty that arises in the case of quasilinear or fully nonlinear 
equations like (|1.5p is that the only well-posedness result which is known in the vicinity of the 
solitary wave is a local well-posedness result in H'^ for s sufficiently large which comes from the 
high order energy method. Note that the global or almost global existence results like the ones 
of [291 So] ) [m [I2] are obtained in a regime where solitary waves do not exist and that the local 
existence results for rough data as obtained in [2] still require a regularity much higher than the 
one which is controlled by the Hamiltonian. 

This makes the perturbative analysis more delicate and requires new ideas with respect to the 
semi-linear setting. 

We thus consider two solitons Qd {x — cit) and Qc2 {x — h — C2t) of (jl.Sp with ci < C2. We suppose 
that ci and C2 satisfy (II. 7p with suitable choices of the small parameters £1^2- We also suppose that 
/i > is large enough. We define 

M(t, x) := Qci {x -cit) + Qc2 (x -h- C2t) (1.8) 

as the two-soliton function. We will focus on the case where each solitary wave is stable in the 
following sense. Under our assumptions (jl.7p on the speed c of a solitary wave, it was proven in [21] 
that for sufficiently small corresponding parameter e, the solitary wave Qc is stable in the sense that 
the second derivative of the Hamiltonian at the solitary wave restricted to a natural co-dimension 
2 subspace is positive. We shall assume that the speeds ci , C2 are such that this property is verified 
(see Proposition 13. 6p . Our main result reads: 

Theorem 1.2. Let us fix s > 0. Suppose that the speeds ci < C2 satisfy ()1.7p with parameters £1, 
£2- Define M by (jl.Sp . Then there exists £* such that for £±,£2 G (0,e*] and h sufficiently large, 
we have that there exists a (semi) global solution U{t) = (rj, c^)* to the water-wave system lll.5\) 
satisfying 

U-M eCb{[0, 00); H'{R) xH'{R)) 

and 

lim \\U{t) - M{t)\\H- =0. 

The assumption that £1^2 are sufficiently small will be only used in order to ensure the exis- 
tence of smooth exponentially decreasing solitary waves (given by Theorem 11.11 for example) and 
that Proposition 13.61 below holds true. Indeed, for the construction of multi-solitary waves that 
we shall perform, the main ingredients that are required are besides the existence of smooth lo- 
calized solitary waves, the stability property given by Proposition 13.61 of each solitary wave and 
the local well-posedness in i7^ for s sufficiently large of the nonlinear system with the existence 
of a quasilinear hyperbolic system type energy estimate. Instead of the smallness assumption on 
£ in Theorem 11.21 we could assume the existence of the solitary waves (solitary waves can also 
be obtained by variational methods for example, [6]) and that each of them satisfies the stability 
property of Proposition 13.61 

We have focused on water waves with surface tension, nevertheless, since the existence of solitary 
waves is also known (see [13] for example) and since some of them are linearly stable, [23] (note 
that nevertheless an estimate like the one of Proposition 13.61 does not hold in this case, the number 
of negative directions of L would be infinite), it could be possible to perform a related construction 
for water waves without surface tension. 

Finally, let us point out that the assumption that h is sufficiently large is only used in order to 
get a solution on [0, +00 [, an equivalent statement would be to take h = and to get a multi-soliton 
solution on the interval [Tq, +oo[ with Tq sufficiently large. 
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There are two main steps in the proof of Theorem II. 2i The first step is to construct a smooth 
approximate sohition of (jl.Sp that tends exponentiahy fast to M as t goes to infinity. This approx- 
imate solution V^ is under the form: 

N 

U" = M + J2^^ 
1=1 

where each V^ is smooth and verifies the crucial property that it is decaying in i?* like e~''^o('^2-ci)i 
for some eo. Each V solves a linear problem with source term. The existence of some V with 
this decay property will be proven by using the spectral properties of the linearization of the water 
waves system (jl.Sp about each solitary wave. 

Once the approximate solution is sufficiently accurate (i.e for N sufficiently large) which basically 
means that the remainder term in the equation has a sufficiently fast decay in time, we shall 
construct an exact solution as a sum of the approximate solution and remainder term that solves 
a nonlinear equation. The main difficulty is to prove the existence of a solution on [0, +00) for this 
problem having at hand only a local Cauchy Theory in H^ for s large. 

Note that this kind of iterated constructions is related to Grenier's argument [H] in order to 
prove that linear instability implies nonlinear instability for quasilinear equations that has been 
used in [21] . 

The main arguments that we shall use to prove Theorem 1 1 . 2 1 can also be used in order to sharpen 
the transverse instability result proven in [23] and construct for the two-dimensional water-waves 
system that is to say when the fluid domain is 

nt = {{X,z) eR^, -H <z< v{t,x)}, 

a solution on [0, +00) of the system which is different from the solitary wave (and all its translates) 
and converge to the solitary wave as time goes to infinity. The result that we shall prove is the 
following. 

Theorem 1.3. Let us fix s > 0. Suppose that c satisfies (|1.7p . For e sufficiently small there 
exists a global solution U of the 2-D water waves system with initial data Uq satisfying U — Qc & 
Cfe([0,oo);i?''(M2) X H'{R^)). Moreover, one has 

dyUo / (1.9) 

and 

lim \\U{t,x,y)-Qc{x-ct)\\H^ =0. 

i— s>+oo 

We shall recall the formulation of the 2D water-waves system in Section [2l 

Remark 1.4. By the remark after |25| Theorem 1.5] we know that this theorem implies the trans- 
verse instability of the solitary wave. 

This result can be compared to classical results about the existence of strongly stable manifolds 
for ordinary differential equations or semilinear partial differential equations. Results as in Theo- 
rem ll.3l were in particular, obtained for semilinear partial differential equations in [lOl IT] for example 
and also in [25] for the KP-I equation. As previously, the main difficulty for the proof of this result 
comes from the fact that the water-waves system is not semilinear. The proof of Theorem 11.31 also 
relies on the construction of a well-chosen approximate solution and of a remainder that solves a 
nonlinear system. The construction of the approximate solution relies on spectral information and 
semi-group estimates that are contained in [2l]. Consequently, we shall first present the proof of 
Theorem 1 1.3 1 in Section [21 this allows us to essentially focus on the construction of a remainder that 
is defined on the whole time interval [0, +00 [. These arguments will be also useful for the proof of 
Theorem 11.21 



The paper is organized as follows. The next section will be devoted to the various steps of 
the proof of Theorem 11.31 In Section [3l we collect some useful bounds for the Dirichlet-Neumann 
map and we prove the key coercivity property of Proposition 13.61 We shall study in Section [J] 
the error that produces in the equation the sum M of two solitary waves. Then in Section [5l we 
shall construct a suitable approximate solution. This will be the most difficult part in the proof, 
the crucial step will be to prove that the fundamental solution of the linearized equation about M 
has a sufficiently small exponential growth. The final step of the proof of Theorem 11.21 is given in 
Section [6l 

2. Proof of Theorem 11.31 

In this section we consider the two-dimensional water-waves system. The fluid domain at time 
t is defined by 

Qj = {(X, z) G ]R3 I _H <z < r]{t, X)} 

where X = {x,y), H > is a, constant and T]{t, X) is the free surface at time t. We use the Zakharov 
formulation recalled in the introduction [3H117] to write the system for the unknowns i]{t, X) which 
is the free-surface and ip{t, X) which is the trace on the free surface of the velocity potential as: 

Z = ^-P.^P + ^ (GM. + Vx, V.,)' _ ^ ^^^ , V^^ (2.1) 

2 2 l + |Vxr/|^ v^l + |Vxr/|2 

where again b is the coefficient of surface tension, g is the gravity coefficient and G is the Dirichlet- 
Neumann operator. 

Since we study here a single solitary wave with speed c, we change frame (x, y, z) to {x — ct, y, z). 
This leads to a new system 

dtf] = cdxT] + G[r]]ip, 

2 2 1 + \Vxm ^/l + WxW 

We also perform the following change of variables 

r?(t, X) = Hi%j^t, ^X), ^{t, X, z) = cH^i^t, ^X) (2.3) 

and simply note (r/, ip) again as (r/, 99) to have the dimensionless water-waves system 

dtf] = dxTj + G[r]\Lp, 

^ ^ 1|T7 12 , l {G[i]]^ + VxV-Vxrif ^.„ VxV (2.4) 

with 

gH „ b 

Observe that there is a slight abuse of notation, since in (j2.4|) the map G[r]] is defined as above but 
with H = 1. If we note U = (77, (/?)*, the system (|2.4|) can be rewritten as 

dtU = F{U) (2.5) 

with 

dxTj + G[ri\ip 

2 2 i + |Vxr/M vi + Wrnl' 
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The solitary wave Qc of the original system (j2.ip becomes a stationary wave solution for the water- 
waves system (j2.4p that we shall denote in this section by Qsix) = {r]^{x),ip£{x)y or simply as 
Q{x). We want to show that there exists a global solution of system (j2.4p near the solitary wave. 



2.1. The linearized operator. As in [2l], we linearize the water-waves system (j2.4p around the 
solitary wave Q^ = {rji^^ip^y. Let us set 

1 + \dxVer 

and 



PeTj := /3Vx 



(l + (5..r?,)2)^ (l + (a,^^)2)| 



Note that since G[77e](/?e = —(9^%, we obtain that Z^ and u^ are decaying exponentially together with 
all their derivatives thanks to Theorem ll.il {(p^ occurs with a derivative in Ve). The linearization 
of (j2.4p about Qs reads 

dtU = JAU, (2.6) 

where U = {tj, (/?)*, J = [ -. ^ I is a skew-symmetric matrix and 

''"V -dA{ve-iy)-G[r^,]{Z,.) G[r^s] 

is a symmetric operator on L^ x L^. The formula for the differential of the Dirichlet-Neumann 
operator with respect to r] which allows to obtain the above expression is recalled in Proposition 
13.11 (5) below. As in [TTj, we can get a simpler form of the linearized system by the change of 
unknowns 

Vi = 7?, V2 = ip- Z.rj. (2.7) 

We get for V = (Vi, V2)* the system 

dtV = JLV, (2.8) 

where L is still a symmetric defined by 

. _ f-Pe + a + ivs- l)d,Z, {ve - l)d, 

"^ " V -UiVe - !)•) Give] 

Notice that since Qe does not depend on y, the study of system (12. Sp can be simplified by using 
Fourier transform in y. In fact, if for some /c G M, 

V{x,y) = e'^yW{x), 

then 

LV = e'^yL{k)W 
with a symmetric operator L[k) defined by 

Here 

Pe,kU = /3{a4(l + {d,,lle?)--2d^u] - k'il + [d^l^efy'^u], 

and Gf^^k is such that 

G[rje]{f{x)e'^'^) = e'^yGe,k{f{x)). 



2.2. The eigenvalue problem JL{k)U = all . We shall need some results about the spectrum 
and the eigenvalues of JL{k) seen as an unbounded operator on L^(M) x L^(M) with domain 
i?^(M) X H^{M.) which are essentially contained in 



Lemma 2.1. We have the following spectral properties of the operators JL{k). 

• There exists e* such that for every e G (0,e*], the solitary wave Qs is spectrally stable 
against one- dimensional perturbation: a{JL{0)) C iM where a{-) denotes the spectrum. 

• For any k ^ 0, the essential spectrum of JL{k) is such that aess{JL{k)) C iM and JL{k) 
has at most one simple eigenvalue of positive real part. 

• If a is an eigenvalue of JL{k) then so is —a. 

• If a is an eigenvalue of JL[k), k ^0, with non-zero real part then o" G M. 

Note that the combination of above properties of the eigenvalues also yield that there is at most 
one eigenvalue of negative real part for JL{k), k ^ 0. 

Proof. These statements are already contained in Lemma 5.1 in |24j . Note that the first statement in 
the above lemma was obtained as a consequence of the work of Mielke |21j . The only additional point 
is to notice that thanks to the reversibility symmetry of the water waves system and the symmetry of 
the solitary wave one has that if U{x) = (r/(x), (p{x)) is an eigenfunction of JL{k) corresponding to 
an eigenvalue a then U{x) = {r]{—x), —ip{—x)Y is an eigenfunction corresponding to the eigenvalue 
—a (this symmetry of the spectrum could also be obtained by using the Hamiltonian structure). 
The last point is a consequence of the symmetry of the spectrum and the fact that there is at most 
one simple eigenvalue of positive real part. D 

Thanks to Lemma 12.11 and in particular the symmetry of the spectrum, we have the following 
counterparts of |24^ Proposition 5.2] and [241 Theorem 5.3] respectively. 

Lemma 2.2. Consider the eigenvalue problem 

JL{k)U = aU, U £ H^(R)x H^{R). (2.10) 

Then : 

(1) 3K > such that if \k\ > K , there is no solution of (j2.10p satisfying cr < 0. 

(2) 3M > such that if\k\ < K , there is no solution of (j2.10p satisfying a < —M. 

From now on, we shall only consider solitary waves Qe with e G (0, e*] as stated in Theorem 11.31 
so that the above spectral properties are matched. 

Proposition 2.3. For e G (0,e*], there exists a > 0, k y^ and a non-trivial U G H^ x H^ such 
that 

JL{k)U = -aU. 

Finally, we also have as a consequence of the above results: 

Lemma 2.4. For every {ao,ko) satisfying k^ ^ 0, Re a^ < and o"o G a{JL{kQ)), the set 
{{a,k)\a G (t(JL(A;))} near {aojko) is the graph of an analytic curve k i— t- (t(/c) with cr[k) a 
real simple eigenvalue of JL{k). Moreover, we can select an eigenvector V {k) of JL{k) depending 
analytically on k. 

Proof. Suppose that (Tq G cr(JL(/co)) with /cq 7^ 0, Re ctq < 0. From Lemma [2?T| (Tq is necessarily 
a simple real eigenvalue. Let Vq S H^ x H^ an eigenvector, JL(A;o)Vo = ctqVq. Consider the 
map F{V,k,a) = JL{k)V - aV . Then Z)y,<^F(Vo, /co,o-o)(f^,M) = JL{ko)U - a^U - iiVq. Thanks 
to Lemma I2.H JL{kQ) — ctq is Fredholm with index zero and its kernel is one-dimensional. This 
implies that the kernel of Dy^(^F{Vo,ko,crQ) is also one dimensional. Indeed, if {U,fi) is such that 
Dv,aF{Uo,ko,aQ){U, fi) = then {JL{kQ) — (Jq)^U = which, thanks to Lemma \2A\ implies that 
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U = XVq, a G M. This in turn implies that /.t = and thus the kernel of D\/^o--F(Vo, A;o,o"o) is 
spanned by (Vo,0). Once again, thanks to Lemma |2. 11 we obtain that Vq is not in the image of 
JL{ko) — ao (otherwise ctq would not be a simple eigenvalue of JL{ko)). Therefore Dv^aFiVo, ko, o"o) 
is surjective. We are in position to apply the simplest form of the Lyapounov- Schmidt method. 
This ends the proof of Lemma 12. 4[ D 

2.3. Construction of the approximate solution. Let us define for U = (77,93)*, the norms 

\\umEs= Yl \\drdP.d^yU{t,.)h2^^2). 

0<a+/3+7<s 

Proposition 2.5. There exists U^{t,x,y) G r\s>oE^ satisfying 

dtU^ = JAC/° (2.11) 

such that for every s > 0, every eg > one has 

||f/°(t)lb^<c,,eoe-('^«-^°)*, t>0 

where —ao is the smallest possible eigenvalue of JL[k) (and thus ctq is the largest possible eigenvalue 
ofJL{k)). 

Note that we use the notations eq and eo for different parameters. 

Proof of Proposition [275[ The construction is close to the one of Proposition 6.1 in |24j . The situ- 
ation is even simpler here since we need less precise information on the asymptotic behavior. We 
already know from (j2.7p that it is equivalent to solve (|2.1ip for U^ and 

dtV'^ = JLV° (2.12) 

with 

U^ = PV', P=Q J 

for yo. 

First of all, one should locate some negative eigenvalue of JL{k). We already know from Propo- 
sition [2]3] that there exists k ^ such that JL{k) has an eigenvalue —6 < 0. Moreover, thanks to 
Lemma |2. 21 we know that the negative eigenvalues —a of JL{k) can only be found when \k\ < K 
and a < M. 

We try to find out the largest a such that —a is still a negative eigenvalue of JL{k). Define the 
setO = {fe|3 — a ^ a{JL{k)), —a < —6/2}. One can see that O is a bounded non-empty open 
set and that k 1— )• —cr{k) is continuous in Cl. We fix ko, ctq such that 

-ao = -a{ko) = M{-a{k)\ k £ Q} < -5/2 < 0. 

Let us fix eo > and an interval Iq C i^ small enough with /cq G Iq so that — ctq < —cr{k) < —ao + €o 
in Iq by the continuity of A; 1— )■ —cr{k). Thanks to Lemma |2.4| we can can choose an eigenvector 
V{k) depending analytically on k on Iq. By elliptic regularity, we have that V{k) £ II°°. Finally 
since one has a{k) = cj{—k) and V{k) = V{—k), let us set / = Iq U —Iq and 

y°(t, x,y)= I e-"^^^^e'^yV{k)dk, t > 0. 

Then V^ is real and is a solution of (j2.12p . We have for any s, a G N that 



'^°(*>-)ll^.(M2) = c/"e-2<^(^)* ^ \a{k)\''^k'^-\d^,^V{k)\l,^^^dk. 



and hence, there exist numbers Cs^a,to such that for any t > 



This yields similar estimates for [/" = PV^. This ends the proof of Proposition 12.51 D 

Proposition 2.6. For any M > 0, there exists 

M+l 

f/« = [/O + ^ JJ'C/i, W G C°°(M+, F°°(R2)), S€R 

i=i 
such that for every j, one has U^{0) = and the estimates 

\\UHt)\\Es < C,je-(j'+i)('"°-^o)*, Vi > 

for eo sufficiently small (eo < ao/2) for some numbers Cgj independent of t and do the eigenvalue 
chosen in Proposition \2.5[ Moreover, define V^ = Q + dU"" . Then V"" is an approximate solution 
of /i2. 5\) in the sense that 

dtV - TiV) = R^P 
where R"'^ satisfying the estimate 

Proof. We follow the idea of proof of |241 Proposition 6.3]. The Taylor expansion of J- is 

M+2 ^f. 

J^{Q + 5U) = F{Q) + Y, -^D^HQW, ...,[/) + 5^'+''Rm,s{U). 

k=l 
Plugging the expansion of U"" into ()2.5p gives the equations for j > 1 that 

dtW -JAW = Y, Yl -DPT[Q]{U^\...,U^''). (2.13) 

p=2 o<;i,...,ip<i-i 

We note the right-hand side of (j2.13p as S^ . Applying Fourier transform in y to (j2.13p , we get that 

dtW - JA{k)W = S^ , j>l. 
So we need to consider first the equation 

dtU - Jk{k)U = F. 
In order to estimate the solution, we introduce 

\U{t)\ls := Y. (\dtdiUiit, Ol^idR) + l^rsf [/2(t, •)!'. 1 

0<o+/3<s k ^ 



with 



\^\ 



2 

4i 



— rf 



l + \DJ2 



+ \k\^\v\l2( 



LH 



We have the following statement. 



Proposition 2.7. Fix 7 > ctq > 0. Assume that F{t,x,k) satisfies uniformly for \k\ < K the 
estimate 

Y \\drd^,F{t,;k)\\L2<Mse-"'\ t>0. (2.14) 

a+[S<s 

Then we can find a solution U of 

dtU - JA{k)U = F, 
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defined for t > such that there exists constant Cg depending on Mg+so for some fixed sq > such 
that uniformly for \k\ < K , we have 

\U{t,-)\L2 + \U{t,-)\x^^<Cse-''\ t>0. 

Proof. Let V = P^^U, then the equation for U is equivalent to the equation 

dtV = JL{k)V + P^^F. (2.15) 

Let a > be such that 7 > (t > (Tq. By a simple lifting argument, we get from |241 Proposition 6.4] 
and the symmetry of the spectrum pointed out in Lemma |2. II that the semi-group estimate 

holds for some fixed derivative loss sq (we could avoid it, but this is not important in our construc- 
tion). By using again the reversibility of the system, we actually obtain that 

l|e*-^^('=Vo||i2nx,^ < Cse^\'\{\Vo\L^ + |l^o|^^+^o), Vt G E. (2.16) 

Let us choose the solution of (|2.15|) given by 

/oo 

Then, thanks to (j2.16p . we have 

/oo 

with Cs having the claimed structure. This ends the proof of Propoposition 12.71 D 

End of the proof of Proposition 12.61 We can finish the proof by induction on j as in the proof 
of Proposition 6.3 in [23] . Assume that we already built (C/')kj-i whose Fourier transforms with 
respect to the y variable are compactly supported in k and that satisfy uniformly for |A;| < Ki the 
estimates 

\U\t, k)\Es < C,,/e-('+^)(^''-^'')*, l<j-l 
where the | • | e^ norm for functions of t and x is naturally defined as 

\a\<s 

To construct U^ , we first observe that thanks to the induction assumption we obtain as in the proof 
of Proposition 6.3 in |24] (hence in particular by using Proposition 3.9 in [23] to handle the terms 
coming from the Dirichlet-Neumann operator) that 

uniformly for \k\ < Kj. Then we define Uj by 

/oo 

Thanks to the semigroup estimate ()2.16p . we get since 7 = (j + l)(o"o — eo) > ctq > that 

\\W{t,k)\\Es < C,je-(^'+i)("°-^o)* 

and the estimate for U^ follows by using the Bessel identity and the fact that U^{t, k) is compactly 
supported in k. 
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2.4. Proof of Theorem 11.31 Recall that we already defined V"' = Q + 6U"' in Proposition 12.61 as 
an approximate solution U of the water-waves system (j2.5p . In order to get a true solution of (j2.5p , 
we still need to consider a remainder U^ such that 

U = ¥" + U^, 

becomes an exact solution. The system for U is 

dtU^ = TiV + U^) - TiV) - R'^P, t > 
[/^(O) to be fixed later. 

Before solving the system for U^ we need to introduce more notations. For a = {ao,ai,a2) G N'^ 
we note 9° = dt'^d^^d^K For U{t) = {Ui{t), U2(t)f and A: G N we define X^ by 



B(n\ +^ u^ P.^^r.A l„+„v, K'^-'^'l 



\a\<k 
and we note 



H^' 



WWx'^r, = sup \\U{T)\\xk. 
*'^ t<T<T 



We define W'' and WJ^j, by 



\u{t)\\w'' = E Wu{t)\\L,^i^2), ||u||^yfc = sup \\u{t) 

\\^U t<T<T 

|a|<fe — — 



IW 



We will use an approximate sequence of solutions {U""} for water-waves problem to prove that 
there exists a global-in-time solution U (t) of ()2.17p . Let {r„} be a strictly increasing sequence 
such that Tn > and T„ — )■ -|-oo as n — )■ oo. First of all, we define U"'{t) as the solution of the 
water-waves system 



'-ni 



c/"(r„) = 0. ^^--^^^ 

Note that we solve the problem backward in time. For the water-waves system there is no problem 
to do so because of the reversibility. We have local well-posedness for (j2.18p as in [23] (see also [2], 
[5]) |26]). The first step is thus to prove that C/"" is defined on the whole interval [OjT""] and that 
it satisfies an appropriate estimate. This will be a consequence of the following a priori estimate 

Proposition 2.8. Let U"{t) be a smooth solution of i2.18\} on [T, r„] satisfying 1 — ||??"(i)||L°° — 
||?7"(i)||L°° > Co > for t G [T, r„]. Then for m>2, s > 5 and t G [T, T^] we have the estimate 

1 1 Tfi t , Tn t , Tfi 

f'Tn 

with uj : M+ — t- [1, +oo] is a continuous increasing function. 

Proposition 12.81 can be directly obtained from j24j . Theorem 7.1. Indeed, let us define the 
isometry ~ by 

U{t, X, y) = {Ui{T^ - T, -i, -y), -U2{T^ - r, -i, -y))' . 

Then, we note that C/" solves (I2.18P if and only if U'^{T,x,y) solves 

a^[/" = ^(T7«(r) + C7"(r)) - F{V''{t)) + R^Pir), r G [0, T"] 

with the initial data C/"(0) = 0. 

We can now convert the a priori estimate of Proposition 12.81 into the following existence result. 
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Proposition 2.9. Let m > 2. Then there exists M large enough in the definition of V"" and 6 
sufficiently small such that the following holds true. For every T„ > there exists a unique solution 
of i2.18\) defined on [0, T'"] and satisfying 

1 - WvIl^ - ll^"(i)llL- > 0, ||C/"(t)|U™+3 < CM,^<5^^+=^e-(^-^+3)(.o-.o)t^ ^ ^ [o,r„]. 

Proof. From |24t Section 8] for example and the above reversibility argument, we know the local- 
in-time existence for the solution C/"(t) of system ()2.18p going backwards from time T^. We 
want to prove that this solution exists on the whole time interval [0,r„]. By Proposition 12.61 and 
Proposition 12.81 we have that, at least for t close to T„ 

t,Tn \ J^ 

+ CM,™52(^^+3)e-2(*^+3)(-o-eo)t^ . (2.19) 

Define 

T* = inf{r G [0,r„] : Vt G [T, T„], ||C/"(t)||x-+3 < 1, 1 - Wti^l^ - ||r/"(i)||Loo > cq > 0}. 
We deduce from (f2l9]) that 

for t G [T*,Tn]. Now we shall fix the value of Af and impose a smallness condition on 5. Let M 
and 5 be such that 

2(M + 3){ao - eo) > uj{C + CM,m6) 
Such a choice is possible thanks to the continuity of oj. Indeed we can first take M large enough so 
that 2(M -I- 3) ((To — eo) > (^{C) + 2 and then 5 small enough so that \oj{C + CM,mS) — ^{C)\ < 1- 
Therefore, we arrive at the bound 

d - '■^" 

Jt 
Integrating on both sides with respect to time from t to T^ leads to 



,(^_^.(C+CM„^S)t I "||f^n(^)||^^^^^^^ ^^^^^^^2{M+3)g.(C+C7M,™5)t-2(M+3)(<xo-.o)t_ 



Jt 

which gives 

||f/"(i)llx'"+3 <^M,™<j2(*-^+3)e-2(A-^+3)(-o-eo)t (2.20) 

for any t G [T*,Tn]. Now we let 6 small enough such that CM,mS'^^^^~^^' < 1 and 1 — ||?/e||Loo — 
CM,m^ > cq. Then by definition, we have T* < 0, hence the solution t/"(t) for (j2.18p can be 
extended to the whole time interval [0, T„] with the claimed estimates. D 

We can now finish the proof of Theorem 11.31 bv invoking some standard compactness arguments. 



Step 1. Existence of the global solution U{t). Thanks to Proposition 12.81 and Proposition 12. 9^ 
we get the approximation sequence {[/"} of solutions of (j2.18p which satisfy 

||f/"(t)||x-+3 < CM,„<5^^+3e-(^^+3)(-o-^o)t^ Vt G [0,r„]. 

Let ip G C^(-l/2, 1/2) be a bump function such that V(a^) = 1 for x G (-1/4, 1/4). We extend 
U^{t) as zero for t > T^ and we set 

U^{t) = i;{t/Tn)U^{t), t > 0. 

12 



Then, we have 






and thus new sequence {W^it)} satisfies 

and 

||9iC7"(t)||^„+3xH™+5/2 < CM,,n5*^+'e-(^-^+3)(-o-^o)i^ Vt > 0. 

By a standard compactness argument, we obtain that there exists a subsequence {U"'''{t)} and 

U^{t) G L°^([0,oo), //"+^ X i/™+V2) 

such that 

^-fe^[/« in CioM+,Hl^+^xH'^^+'/^) as n^ ^ oo 

and 

r''WllH"^+4,<^™+7/2 < CM,m<5*'+=^e-(^^+3)(<xo-.o)t vt G [0,oo). 

Moreover U^{t) is the solution of plTl) since ?/"'=(*) is a solution of (f2Tfll on (0,r„j4). Going 
back to water-waves system ()2.5p we deduce that there is a global solution U(t) = V^it) + U^{t) 
for system (j2.5p . 



Step 2. Behavior when t — t- +oo. By the definition of V"' we have 

Af+l 

C/(i) = y"(t) + c/«(t) = Q + J^ 5^+i[/^(t) + U'^it), 

j=0 
and from Proposition 12.61 and (J2.2U|) we know that for any t £ [0, oo) 

with s > 2. This yields 

WU^mH^ < CA/,.e-(j'+i)(-o-^o)*, ||C/^(t)||H» < CA/,.e-(*^+3)(-o-eo)t^ Vi G [0,oo) 

which shows that 

lim \\U{t)-Q{x)\\H^ =0. 

Step 3. It remains to check the condition (jl.9p for [/. Since 
we have at i = that 

M+l 

[/(O) = Q(x) + (5 ^ 5^C/J(0) + C/^(0) 

j=0 

and 

Af+l 

9,([/(o)) = sdyiu'm + Y^ 6^^'dy{wm + dyiu^'m. 

i=i 
We know from the definition of U^ that dy{U^{0)) / 0. We can use Proposition 12.61 and step 1 to 
get that 

This yields that dy(U{0)) ^ for 5 possibly smaller. This ends the proof of Theorem 11.31 D 
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3. Preliminaries for the proof of Theorem 11.21 

3.1. The Dirichlet-Neumann operator. Let us recall that the Dirichlet-Neumann operator G[r}] 
is defined by 

where n is the unit outward normal vector and ip solves 

Aif = 0, in — H < z < rj{x) 

(p\z=r, = 1p, dn^p\z=~H = 

We can rewrite the elliptic problem on —H < z < rj{x) as an equivalent problem on a flat strip 

S = Rx [-1, 0]: 

Va;,2 • P'^x,z<f = 0, in 5 ,g ^. 

(f\z=0 = ip, dz(p\z=-i =0 
where (p{x, z) = (p{x, Hz + [z + I)/?) and 

With the notation d^ = (0, 1)* • PVx,z, one can see that the D-N operator associated to the above 
problem can be written as 

G[ri]ilj = d^Lp\z=o = -dxVdx'f\z=o H ftv dz(p\z=o- (3.2) 

Let us recall the following properties: 

Proposition 3.1. For rj G H°^{M.) with H + 7] > cq > 0, we have 

(1) G[rj] is symmetric on L'^{M.): 

{G[ri]u,v) = {u,G[r]],v), Vn,fGF^(IR) 

(2) There exists c > 0, C > such that for every u G H^{W) 

\{G[rf^u,v)\ < C|*Pn|i2|«Pz;|i2,V'u, v G H^{R) (3.3) 

{G[r]]u,u) > c\^u\l2, Vu G H^{R) (3.4) 

where *P is the Fourier multiplier 

(3) the linear operator G[r]] : //^"'"^(M) — t- H^{M.) is continuous for every s G M 

(4) We have the following commutator estimates 

m,G[v]]u\^i < Cs^iulHs, Vug //^+5(M), 

\{fdxU,G[v]u)\ <Cf\^u\L2, VnGif^(M),V/G//°°(E). 

(5) We have the explicit expression for the derivative of the Dirichlet-Neumann operator with 
respect to rj: 

R,{G[i^]u) ■ C = -GMiCZ) - dxivC) 

with 

Z = Z[r], u] = — ,„ '"i^'' , V = v[r], u] = d^u - Zd^rj. 

J- + \C'x'n\ 
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(6) Finally, for s > 1/2, we have that 

j 

■j=i 

For the proof of these statements we refer to [T7] or [23] section 3. We have assumed that r] is 
smooth since we shall use this proposition when 77 is a soHtary wave. Most of the estimates are 
actually still true when rj only has limited Sobolev regularity. 

In order to perform our construction of multisolitons, we shall also need results about the action 
of the Dirichlet-Neumann operator on localized functions. 

Proposition 3.2. Assume that ip E C^(M) has an exponential decay: 

3d>0, VaGN,a> l,3Cay x eR,\d^ip{x)\ < C^.e"'^!^!. (3.5) 

Then for rj G H°^{M) with H + rj > cq > 0, G[ry]'0 also has an exponential decay, that is, for any 
a € N, there exist a constant Ca depending on a and < e < d independent of a such that for 
every x G M, 



Remark 3.3. We only assume that the derivatives of the function ij) are exponential decaying while 
the function itself is only hounded as it is the case for the solitary waves (see ()4.2p below). Note 
that this still yields that G[7]]^l' is exponentially decaying. The heuristic reason is that the Dirichlet- 
Neumann operator behaves like a derivative. Again, in Proposition 1 3. S\ we are not interested in the 
way the estimates depend on the regularity of the surface since we will use it for solitary waves. 

The result of Proposition 13*.^ uses in an essential way that the fluid domain is bounded in the z 
direction (we use a Poincare inequality). The statement of Proposition [3^ does not hold in the case 
of an infinite bottom because of a singularity at the low frequencies which affects the propagation of 
the exponential decay. 

Proof of Proposition \3.SX We use as in [23] the decomposition 

/^{x^z) = ipQ{x,z) +u{x,z) (3.6) 

with ifQ that solves the elliptic problem 

-A^.,^C(5o = 0, {x,z) G 5, fo{x,G) =■0, dz^o{x,-l) = 0. 



This allows to transform the elliptic problem (j3.ip for (f into an elliptic problem with homogeneous 
boundary conditions for u: 

u\z=o = 0, dzu\z=-i = 

At first, we shall study the decay properties of ipQ. We first observe that fo is given by the explicit 
formula 

_, , cosh(^(z + l)) ^^ ^ , ^ 

cosh(^j 

where ^ stands for the Fourier transform in the x variable. From this expression, we get in particular 
that 

Md.v>omz) = ^^^^^%^-F.(9.^)(6. 

cosh(4j 
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The exponential decay will follow from Paley- Wiener type arguments. Since dxip and its derivatives 
have exponential decay, we get that J-x{dxtp) has an holomorphic extension to |Im,^| < d and by 
integration by parts that it satisfies for every 5 G (0, d) the estimate 

\j'x{dx^Gmz)\ < Y^|w' ^^' ^' Iim ^\<S,ze [-1,0] 

for every A^ G N. Since S, i— )• — ^^^^^/g^ — - has an holomorphic bounded (uniformly in z) extension 
to |Im 1^1 < 5 for any 6 E (0, 7r/2), we can use contour deformation to write that 

Jr cosh(^) 7/^5=5 cosh(0 

for any 6 such that \6\ < Min (d, 7r/2). This yields by choosing (^ = 2esign x, with e sufficiently 
small the estimate 

In a similar way, we get from (j3.8p that 

sinh(f(z + l)) -._ sinh (£(z + 1)) 1 , 

*'^»<^-'> = cosh,; ^ "-'fl = coshf '^ -■^.a.^OK) 

and hence, since dxip and its derivatives have exponential decay, the same arguments as above yield 
the estimate 

The estimates for higher order derivatives can also be obtained from the same arguments, we find 
in the end: 

V/3, |/3|>1, \d^^,Mx,z)\<Cpe-^^\^\, y{x,z)eS. (3.9) 

It remains to estimate the solution u of ()3.7p . Let us define v{x, z) such that 

u(x, z) = e~'"^'^^v{x, z), (x) := (1 + x^) 2 

with e > small enough and to be fixed later. One has the elliptic problem for v{x, z) 

-V.,z ■ (.Py.,zv) - [Vx,z ■ Pyx,z, e^<-)]e-^<-)^ = e'^^)Vx,z ■ {PV.,z^o), m S 
v\z=o = 0, dzv\z=~i = 

We shall first estimate Sobolev norms of ^(x, z). 

1) Lower-order elliptic estimate for v. By integration by parts, we get that 

P^ x,zV ■ V x,zvdxdz 
s 

= I ([V,,, • PV,,„ e^<">]e-^<">?;) vdxdz + / e^<">V,,, • [PV x,zm) v dxdz 

:= A1+A2. 
For the left hand side, we have by the assumption on 77 that 

P'^x,zV ■ Vx,zvdxdz > c\\Vx,zV\\ 
s ' ' 

for some c > independent of e. Here and in the sequel the norm || • || is the L?'{S) norm. Next, 
we can compute 

[V,,.-PV,,., e<^^] = {H + 'n)[dle<^^]-2{z + l)dxr,[dx,e^^^^]dz 

which yields 

l^il <C(e||V^,^7;||+e2||t;||)||t;||. 
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(3.10) 



For the second term in the right hand side, we can use (j3.9p to get 

|^2| <C,\\v\\ <e^\\vf + Ce, 

where the last inequality comes from the Young inequality, for some harmless number C^ (that 
depends on e). Summing all these estimates up one gets that 

\\Vx,zvf < Cie^Wvf + e\\V^,,v\\ \\v\\ + C,) 

To conclude, we note that since v\z=o = and S is bounded in the z direction, we have the Poincare 
inequality: 

Ikll < c'liv,,,^;!!. 

Plugging this into the above estimate yields 

||Va;,2U|| < C, \\v\\ < C. 

by taking e sufficiently small. 

2) Higher-order estimates for v. These estimates will follow from an induction argument and 

standard elliptic regularity theory. Indeed, we can write the elliptic problem ()3.10p under the form 

-V^,^ • (PV^,,t;) = F, (x,z)e5, v(x,0)=0, a,t;(x,-l)=0. 

From standard elliptic regularity theory, we have that for s > 

lbllH-+2(5) < C[\\v\\h1(^s) + \\F\\h^(^s))- 

By using the estimate (13.9P and a standard commutator estimate, we get 

||i^||H«{5) < C{l + ||f^||H«+i(5)). 

Consequently starting from the H^ estimate that we have already proven, we get by induction that 

for every s > 0, ||'y||H''(5) ^ Cg- By Sobolev embedding, we thus obtain that e'^^^'^^ '"^ u and all its 
derivatives are bounded in S that is to say: 

V/?, \dlXx,z)\ < C^e-^l^l, y{x,z) £ S. (3.11) 

To end the proof of Proposition 13.21 it suffices to combine the expression (13. 2p of the Dirichlet- 
Neumann operator (note that it always involves a derivative applied to <p) with the decomposition 
(j3.6p and the estimates ()3.9p . (13. lip . This ends the proof of Proposition 13.21 D 

We shall also use the following corollary. 

Corollary 3.4. Let ip G C^(M) with an exponential decay property as in Proposition \ 3.^ Then 
for rj G if°^(]R) such that H + rj > cq > and every a £ R, G[r/]('0(a; — a)) also has an exponential 
decay, i.e. there exists constant < e < d independent of a and a such that 



\d^GMii;{--a))ix)\<Co^e 



-e\x—a\ 



Proof. Let us introduce the translation operator {Taf){x) := f{x — a). The result follows by ob- 
serving that G[?7](raV') = 7"a(G'[r„a^]V') ^nd by using Proposition 13.21 D 
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3.2. Linear stability properties of the solitary wave. In this section, we study the hneariza- 
tion of the water-waves system about a sohtary wave Qc{x — ct) given by Theorem 1 1.1[ The main 
results of this section (in particular Proposition 13. 6p are essentially contained in [21] . For the sake 
of completeness, we shall give proofs that use a slightly different framework which is more adapted 
to our purpose. 

It is convenient here to go into the moving frame by changing x into x — ct. As in section I2.H 
the linearized equation reads 

dtU = JAcU 

with 

_ f-Pc + g + Z,G[r],]{Z,-) + ZAv, {vc - c)d, - Z,G[ri,] 
'"V -d.{{vc-c)-)-G[r^,]{Z,-) G[7?,] 

and Zc = Z[Qc\, Vc = v[Qc] are defined in (5) Proposition 13. 11 The operator Pc is the second order 
elliptic operator defined by 

Pc = bd,. 



dx- 



[I + {3x^0?)- 
We can also get a simpler form of the linearized system by the change of unknowns 

V = RcU, Rc=( }^^ \ \ (3.12) 

which yields 

dtV = JLcV, (3.13) 

with a symmetric operator Lc defined by 



^ ^ f-Pc + g + {vc- c)dxZc (vc - c)d, 
-dx{{vc-c)-) G[n" 



X 

■c' 



The two operators are related by the property 

L, = {R^yA,R-\ (3.14) 

Note that in this section, we see c as a more natural parameter than e for the solitary wave and 
the objects depending on it since we have not written the system in a non-dimensional form, but 
that Lc is conjugated to the operator L(0) studied previously via the scaling transformation (j2.3p . 

Thanks to Proposition l3.lt the quadratic form associated to Lc is naturally defined on the space 

1 1 

X^ = H^{]R) X H^ (M) where H^ (R) is a modified homogeneous Sobolev space defined by 



1 

-2/ 



H^{M.) = {ue S'{R), ^iu G L^{R)}. 

Remark 3.5. Note that with our definition of ^, we have that H^ (R) C L^^^(R). Indeed if 

u G H^iR), then v = *Pn G L^ . Let us choose x(0 a smooth compactly supported function such 
that X = ^ in the vicinity of zero. Then since 

{l-x{D))u = ^-'{l-x{D))v, (3.15) 

and that ^^) is bounded, we get that (1 — x{D))u £ L'^{R). Next we note that 

XU = C+ r{l-dl)-^{xv)dy (3.16) 



for some constant C. Since (1 — d^)i{xv) G L (R), this yields that x{D)u G Lf^^iR) and hence 
that 

u = xu + {l-x)n<^Llc{R). (3.17) 
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1 
On H^ (M), we shall use the semi-norm |n| . i = |*!pn|/,2 and hence on X^, we set 

\U\xo = |C/iUi + |C/2|.i. 

We could make X^ a Banach space by taking the quotient by the linear space 

0^0 = {(0, A), AG M}. (3.18) 

Nevertheless, we refrain from doing it since it is very convenient for us to have that X^ is a subspace 
ofLL(or50. 

Proposition 3.6. Suppose that c satisfies ()1.7p . There exists e* such that for every e G (0, e*], we 
have for Qc = {rjc, ^c) the corresponding solitary wave of speed c that there exists C > such that 
for every U = {Ui, U2) £ X^ such that {U, JRcd^Qc) = {Ui,d^r]c) = 0, 

{L,U,U)>C-^\U\\o. 

The orthogonality condition that appears in this proposition is with respect to JRcdxQc because 
of the relation (j3.14p since we are dealing with the operator Lc- If we translate the result into a 
positivity property in terms of Ac, we recover the natural condition {JdxQc, U) = 0. 

Remark 3.7. A more natural orthogonality condition related to the general Grillakis-Shatah- 
Strauss framework [l5] would be to show that if 

{U, JRcdxQc) = {U, RcdxQc) = (3.19) 

then 

iLcU,U)>C-^\U\j^o. (3.20) 

Nevertheless, note that the condition (U, RcdxQc) = is not really well defined on X^ because the 
scalar product {U2,dx4'c) would not be uniquely defined in the quotient since ipc has different limits 
at ±00. Using this remark differently, we observe that for U under the form 

U = XiRcdxQc + X2{0,iy, 

for every Ai, one can find some A2 such that {U, JRcdxQc) = {U, RcdxQc) = 0. This yields that we 
can always find some U such that LcU = and that satisfies the orthogonality conditions ()3.19p 
but that U ^ 0^0 since one can always chose Ai 7^ 0. Therefore, (I3.20p under (I3.19P is false. 

Proof. We first prove the following weaker version of the statement. 

Lemma 3.8. Let [/ = (f/i, C/2) G X°, U ^ O^o and such that 

(U, JRcdxQc) = iUudxVc) = 0- 
Then {LcU, U) > 0. 

Proof of Lemma \3.8[ Lemma 13.81 is a variation on \2A\ Proposition 4.8] with different orthogonal- 
ity conditions. It turns out that the orthogonality conditions used in |24l Proposition 4.8] are 
not appropriate for our purpose here. Using the rescaling ()2.3p . we already know thanks to \2A\ 
Proposition 4.8] that: 

Proposition 3.9. There exists a co-dimension two subspace X^ of X^ and a constant cq such that 
for every U G Xf, 

{LcU,U) > co\U\lo ■ 
As a consequence, thanks to (j3.14p and using that R^^X^ C X^ , we obtain that there exists a 
co-dimension two subspace Xg of X^ and a constant cq such that for every U G Xg, 

(AeC/,C/)>Co|C/|^0. 
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We shall now prove Lemma [3^ by using Proposition 13. 9i Note that because of (j3.14p and using 
that [R~^yj = JRc, it is equivalent to prove that for every U ^ 0^0 and such that {U, JdxQc) = 
{UijdxTjc) = 0, one has {AcU,U) > 0. Recall that the solitary wave Qc is a critical point of the 
Hamiltonian H(r],ip), defined by 

Hiv, 'P^ = ir {G[r]]^ip + gri" + 26(^1 + {d^riY - 1) - Icr^d^^ 

Thus 'VH{Qc) = 0. Differentiating the last relation with respect to c and x respectively gives 

AcdcQc = JdxQc, AcdxQc = 0. (3.21) 

By contradiction, let us assume that there exists y = (2/1,2/2) G X^, y ^ Ojfo and satisfying 

(2/, JdxQc) = {yi,dxVc) = 0, (Ac2/, y) < . 

Set Y = spandi, dxQcT dcQc)- We claim that dimy = 3. In order to prove this claim, we can use 
the following lemma. 

Lemma 3.10. Under the assumptions of Proposition [STSl dxQc oind dcQc o-re not co-linear and 

{dcQc,JdxQc)<0. (3.22) 

We will give the proof of Lemma 13.101 later. Let us use it to prove that the vector space Y is 
three-dimensional. If we suppose that 

y = adxQc + ^dcQc, a,/3GlR 
then by taking the scalar product (the distributional duality) with JdxQc and (^^^^/ciO)*, we get 
= (y, JdxQc) = P{dcQc, JdxQc), = (2/i,9^7?c) = a\\dxVc\\'i2 + P{dcr)c,dxVc)- (3.23) 

Using Lemma 13.101 and (I3.23p . we obtain that a = /3 = 0. Therefore dimy = 3. 

Now, a similar argument shows that Y n O^o = {(0, 0)*}. Indeed, it suffices to use that (0, 1)* is 
orthogonal to JdxQc and (S^^oO)* and to take the scalar product with these vectors in a relation 
of type 

A12/ + X2dxQc + X^dcQc = (0, 1)*, Ai, A2, A3 G M. 
Next, for ^J-i, fJ-2, fJ-s £ ^, we can write by invoking Lemma 13.101 and (I3.21|) . 

{Acifiiy + ^i2dxQc + fJ'3dcQc),fJ'iy + fJ'2dxQc + fJ'sdcQc) = /^?(Ac2/, y) + ulidcQc, JdxQc) < . 

Therefore {AcU,U) < for every U £ Y. But since dimy = 3 there exists a nontrivial z £ Y 
such that z £ X2 where X2 is the space involved in the statement of Proposition 13. 9[ Therefore 
\z\xo = 0. This implies that z E O^o. But z £ Y and Y n O^o = {(0,0)*} which implies z = 0. 
Contradiction. This ends the proof of Lemma 13.81 D 

Proof of Lemma \3.10[ Thanks to [3] (see also Theorem II. ip . we have by setting X = x/H 

|(x) = -e\h~'(—^^) + O(.^e-I^l) (3.24) 

H ^2(/3- 1/3)2^ 



^c /„_^ ^_^o 1 /oNi^i,-2/ ^^ \ . /o/-3^ 



-(x) = -2e(/3 - l/3)2th-^ -^ + 0{e'), (3.25) 

cH ^2(/3- 1/3)2 



and 



where /3 = jr-^ and e = y ^ — 1- Then dc£ < and using that 

/oo 
r]c{x)dx(pc{x)dx (3.26) 

-00 

we deduce <^^?27\i by substituting ([HT^ and (pl^^ in (pl^ . Finally, one also deduces from ([^^ 
and (|3.25p that dxQc and dcQc are not co-linear. This completes the proof of Lemma 13.101 D 
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Let us now come back to the proof of Proposition 13.61 We argue by contradiction. Suppose that 
there exists a sequence (JJ") such that 

|C/"|xo = l, hm(L,f/",C/") = 0, {U'',JRcd,Qc) = {U^,d,7],)=0. (3.27) 

n— >oo 

Then (up to the extraction of a subsequence), there exists Ui G H^ such that (C/f ) converges weakly 

_ . 1 

in H^ to f/i, and, by using Remark 13.51 there exists IJ2 € H^ such that U2 converges weakly in 
^foc towards U2 and ^f/^^ converges weakly in L^ towards '^1)2 ■ Next, we set (j = (C/i, C/2) G -'^'^■ 
We have that U satisfies the same orthogonality conditions as C/": 

{U, JRcd^Qc) = {Ui,d^7]c) = . (3.28) 

Indeed, the second assertion is a direct consequence of the weak L^ convergence of C/" to Ui. In 
order to prove the first one, we write 

= ([/", JRcd^Qc) = (f/r,a,9?c - Z.d^r^,) - {U^.d^r,,) . 

From the weak L^ convergence of (C/") to C/i, we obtain that 

([/", dxipc - Zcd^ric) -^ {Ui,dxipc - ZcdxT]^) 
and by observing that 

{U^.d^r,,) = -[W2, (1 - dl)'^%), (3.29) 

since (1 — d^)^'^7]c G L"^, we also get by weak L^ convergence that 

m,d^r^,) ^ -fe[/2,(l - diy/%) = (C72,a,.7?,) 



(the last scalar product is well defined thanks to the exponential decay of dxijc)- This proves that 
U verifies dS^HD- 

Let us now complete the proof of Proposition 13.61 Let us set 

£c = -bdx {Ccdx ■) +g + {vc- c)dxZc, Cc = (1 + \dxVc\'^)~"', 
we shall first conjugate to get a leading order part with constant coefficients: 



hence 

Let us set 
so that 



ly "- vCc 



rricScimc-) (fc - c) rricdx 

dx{{vc - c)mc-) G[r]c\ 






{L.U'', [/") = {L.V'', V) . (3.30) 

Note that we still have that V{^ is bounded in H^ . We can thus assume (after extracting a subse- 
quence) that y" converges strongly in Hf^^ for every s < 1 towards Vi = l/rricUi. For convenience, 
we also set V2 = U2 so that V2" converges weakly in Lf^^ towards V2 and that ^^^2" converges weakly 
in L^ towards ^V2. 

Next, we observe that we can write the decomposition 

i. = L.+A',L.= (-«|+» -«^) (3,31) 
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and 

K = f ^'"''~ '^)^cdxZc + g{mc - 1) + bmcdximcCc)dx + hmcdx[C,cdx'mc ■ ) {vcruc - c{mc - l))dx 
^ V -dx{{vcmc + c{l-mc))-) 

For Ki, which is a relatively compact perturbation, we obtain 

lim (KiV, y") = (KiV, V). (3.32) 

Indeed, we can use the same trick as in (j3.29p . the exponential decay of Vc, dxZc, rUc — 1, dxTUc, 
d^rric and the fact that V{^ converges strongly in fff^^ for every s < 1. 

To analyze Li, we use the following factorization of its associated quadratic form: for every V 

(LiF, V) = {{-bdl + g- ^M-^)Vi, Vi) + {M{dxV2 - cM-^i), S^Fs - cM-Vi) 

where 

M = -d-^G[i^,]d-\ 

Note that M is a well-defined operator (of order —1), positive and invertible thanks to Lemma 4.5 
of [211 . Let us set 

Mo = -5-iG[0]5~\ 
we can rewrite 

(LiF, V) = {LiV, V) - c^{M-^ - M^')Vu Vi). 
with 

{LiV,V) = {{-bdl+g-cH'I^^)Vi,Vi) + {M{dxV2-cM~^V,),dxV2-cM-^V,). (3.33) 

Note that M — Mq is a compact operator on H^ (see the remark at the bottom of p. 295 and the 
top of p. 296 in [21]), consequently, we have 

lim ((M - Mo)V{\ V{') = ((M - Mo)Vi, Fi) . (3.34) 

n 

To pass to the limit in (LiV^jV^), we note that Li is a non-negative operator. Indeed, M is 
nonnegative and the first part of Li is a Fourier multiplier with symbol 

"*' = '«^ + " - =^i^ = f («^«^ + " - i^) 

which is positive since /? > 1/3 and a > 1 thanks to ()1.7p . Consequently, we get from the weak 
convergence properties that 

liminf(Liy'',y'') > {LiV,V). 

Gathering the previous transformations, we thus get that 

liminf (Lc?7'',C/") =limmf (L^V^V) > (LiV,V) - c^((M - Mo)Vi,Vi) + (KiV,V) (3.35) 

= {L,V,V) = {LM,U). (3.36) 

Consequently, thanks to (j3.27p . we find that {^LJJ , f/) < Since JJ verifies (j3.28p . we get thanks to 
Lemma 13.81 that U G Oxo that is to say U = (0, A)* for some constant A S M. Since Ki is such that 
{Ki{{^, A)*), (0, A)*) = 0, we obtain by using again (ISlSOD . (IHISTD . (13:32]) . (l33iD that as n ^ oo, 

(LeC/", U"") = {LiV, V") + o(l). 

Moreover, since the symbol m(^) is positive and the operator M is non-negative we get from ()3.33p 
that for some cq > independent of n, we have 

{L,U'',Un>co\\Vm,+o{l). 
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This yields that V{^ converges strongly to in H^. This is sufficient to also obtain that 

n 

and hence we find by using again (j3.30p . (|3.3ip that 

{L,U^,Un = {G[vc]U^,U^)+o{l). 

Thanks to ()3.4p in Proposition 13. H this implies that ^U2 converges strongly to zero in L^. We 
have thus obtained that lim„ Hf^^Hxo = which contradicts the assumption that Hf/^Hxo = 1 in 

This ends the proof of Proposition 13.61 D 

From Proposition 13.61 one can get that 
Proposition 3.11. For every U € X^ there exists a unique decomposition 

U = aJRAQc + (iRcdxQc + V (3.37) 

with V € X^ such that 

(y, JRcd^Qc) = {Vi,d,r],) = 0. (3.38) 

Moreover, there exists cq > and C > such that for every U € X^ written under the form (|3.37p . 
one has 

(L,C/,C/)>co|y|^o-C|a|2. 

Note that in the decomposition (|3.37p . V is not orthogonal to the RcdxQc- This decomposition 
has better properties than the orthogonal decomposition that one would get from proposition 13.6 
by choosing V orthogonal to JRcdxQc and [dxTjc^Y ■ The reason is that RcdxQc is in the kernel 
of Lc while {dxTjc, 0)* is not. 

Proof of Proposition \3.11\ If 

U = aJ RcdxQc + PRcdxQc + V (3.39) 

with V satisfying (|3.38p then a and /3 are necessarily determined by 

{U, JRcdxQc) (^,(9,.??e,0)*) 

\JRcdxQc\\2' {RcdxQc, (dxTlcOy) 

and thus are well-defined since {RcdxQc, {dxi]c,^Y) = |<9x??c|^2 / 0. This proves the uniqueness 
and one directly verifies that with a, (3 defined by (j3.40p . the function V in (j3.39p satisfies the 
orthogonality conditions (j3.38p . 

Next, since LcRcdxQc = ^cdxQc = 0, we get by using the decomposition p.39p that 

{LcU, U) = a\LcJ RcdxQc, JRcdxQc) + 2a{V, Lc J RcdxQc) + {LcV, V). 

Therefore, by using Proposition 13.61 for V, we obtain 

{LcU,U) > ^\\V\\^j,o - Ca\ 
This completes the proof of Proposition 13.111 D 

As a simple corollary, we can get stability estimates in X^ for the linearized equation 

dtU = JLcU, U/t=o = U^. (3.41) 

Corollary 3.12. There exists C > such that for every U^ G X°, the solution ofU{t) of (|3.4ip 
satisfies the estimate 

\U{t)\xo<C{l + \t\)\U%o. 
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a = I ^„ ^ ^ 12 , /3 = ^D a ^ ^a „ n\t\ ~ ^{JRcdxQc, {dxr]c, 0) ) (3.40) 



Proof. From the explicit expression (j3.40p for a, we get dta = 0, and therefore 

Ht)\ = |a(0)| < C{\\U?\\l2 + \{Uld,Vc)\) < C\U\o . (3.42) 

Since U solves (j3.4ip . we observe that dt{LcU, U) = 0. Therefore, we get from Proposition 13.111 and 
(1021) that 

\V{t)\j,o < C((L,[/0, C/O) + |a(t)|2) < C\U%o . 
Moreover, from the explicit expression (|3.4U|) for /3, we find that 

{JL,U, {d^r,,, 0)^) _ {JL,{aJR,d.,Qc + V), {d^Vc, 0)*) 
*^ {Rcd^Qc,{d..r],,oy) {Rcd-,Qc,{dx^c,OY) 

Hence 

\m\ < C\U'>\xo + C f\\a{s)\ + I V(s)|;,o)ds < C(l + |t|)|C/°|xo. 

^/o 

This completes the proof of Corollary 13.121 D 

4. Error produced by a sum of solitary waves 

We shall now begin the construction that will allow to get Theorem II. 2i We denote the two 
different soliton solutions for the water-wave system as (jl.5p as 

Qci{x - cit) = Qi{x - cit) = {r]i{x -cit),(pi{x-cit)y 

Qc^ix - h- C2t) = Q2{x - h- C2t) = {r]2{x - h- C2t),(p2{x - h- C2t)y, 

and we set 

M{t,x) = Qi{x -cit) + Q2{x- h- C2t) := {r]Mit,x),ipM{t,x)Y 
to be the superposition of the two solitary waves. Sometimes we shall simply write Qi = {i]i, ipiY 
and Q2 = {r]2^'^2Y for convenience. Thanks to Theorem I l.l) we have 

Bd > 0, Va, 3C«, Vx G M, |aj%r/i| < C<,e-'^l^-'='*-(*-^)''l, i = 1, 2 (4.1) 

and 

3d>0, VaGN, a> 1, 3C„, VxGM, \dl^^i\ <Cc,e-'^\^-^^^-^'-^^'''\ i = 1, 2. (4.2) 

Note that ifi is bounded but not exponentially decaying, nervetheless, derivatives of ipi are expo- 
nentially decaying. In this section, we shall establish that M solves the water-waves system (|1.5p 
up to a small exponentially decaying term: 

Proposition 4.1. The two-soliton M{t,x) solves: 

dtVM = G[r]M]ipM + Ri 

dtVM - -2\0x^m\ + 2 i+|a.r,Mr^ 3"^^^ + ^^^ \ ^i+\d.m,? ) + ^2 

where the remainder RM{t,x) := (i?i,i?2)* has an exponential decay in time, that is, there exist 
constants Cg and eo > such that for any s > 

Let us recall the notation 

\U{t)\Es = Y. K.U\l^- 

\a\<s 

Note that in the sequel we use again both e and e for different parameters. 

In order to prove Proposition 14.11 the main difficulty is to study the interaction of the two 
solitary waves via the Dirichlet-Neumann operator, i.e we need to study for example G[rji\ip2, thus 
the crucial ingredient will be Proposition 13.21 
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4.1. Proof of Proposition 14.11 The basic idea is that the interaction between the solitary waves 
is weak because they are localized and far way with different speeds. We shall use many times in 
the proof the following elementary lemma: 

Lemma 4.2. For C2 > ci, e > and eo G (0,e), there exists C > such that for every h > 0, 
t> 0, 



/ 



^-e\x-cit\ ^-e\x-h-C2t\ ^^ ^ Q^-eh^-eo{c2-ci)t 



Proof. It suffices to decompose the integration domain in the three regions x > h + C2t, x < cit 
and cit < X <h + C2t. D 

Let us prove Proposition 14.11 Since Qi and Q2 are two solutions of (ll.Sp . we can sum up the 
first equations of the two systems to get 

dtVM = dtrii + dtri2 = G[7]i]lpi + G[r/2]v?2 

= G[r]^[]LpM + G[t]i]lpi - G[?7m]</2i + G[r]2]'^2 - G[7]M]'f2- 

So we have the first equation for Af : 

dtrjM = G[r]M]^M + Ri 

with i?i = G[rii]Lpi — G[riM]^i + G[ri2\^2 — G[rjM\^2- Next we need to estimate i?i . Using the shape- 
derivative formula for the Dirichlet-Neumann operator (see (5) Proposition 13. ip . we can compute 
that 

G[i]i\ipi - G[r]M]<fi = - Dr,G[rii + sri2\<fi ■ r]2ds 



/o 

[G[r]i + sr]2]{r]2Zis) + d^{rj2{d^Lpi - Zisd^{rji + sr/2))] ds 
'0 

where 

G[r]s]ipi + dxT]sdx(pi 



Z 



1 + \dxris\'^ 
with the notation rjs = i]i + srj2- Next, we can use Corollarv 13.41 to get 

\dr,AG[ris]^i)\<C^e-^\--'^'\. for a G N, 

Indeed, since the solitary waves depend on x — Cit, one can always convert a time derivative into a 
space derivative. We thus get 

la^^^Zi,! < C„e-^l^-"i*l, for a en. 

With these estimates and using Lemma 14.21 we get 

Jo 

< Cse-'°^e-'°^^^-''^^\ for t > 0. 

Similarly we have the estimate for \G[ri2]^2 — G[rjM]^2\E''- Summing these two estimates up leads 
to 

Now we deal with the equation for ipM- From the second equation of system (II. 5p for both Qi and 
Q2 one can write down that 

dt^M = dfipi + dtip2 

= -i^lOxV^Ml +7^ rm^ 12 9r]M + bdx\ = ] + R2 
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where R2 = R21 + -R22 + R23 with 

2 



-R21 = -^{dx^Ml"^ - ^ ■^\9x^i\'^ = dx(pidxip2, 

i=l 



1 {G[r]Af]ipM + OxfudxtlMf' , Y^ 1 (G[?7i]99i + dx(pidxr]if 
^22 - -o 1 , |;^..„„|2 + 2^ 



2 l + |5,r/Afp ^2 l + |a,r/,|2 



i=l 






With the same arguments that we have used for the estimate of R\ , we get that R2 also satisfies 
the same exponential-decay estimate 

\R2\E' < C,e-'o''e-"''(^2-ci)t^ fQ^ ^ > Q_ 

This ends the proof of Proposition 14.11 D 

5. Construction of an approximate solution 

If we take 6 = e~^"^ > 0, 6 will be small enough if we choose h > large enough later. The 
remainder Rm in the system for M{t, x) can be rewritten as Rm = SRm with 

\Rm\h- < C.e-^o^'^^-^i)*, for t > 0. (5.1) 

Let 



V{t,x) = Y,S'Vi{t, 



N 

-. :X) 
1=1 

with unknowns Vi{t, x) to be constructed later. We want to show that U"'{t, x) = M{t, x) + V{t, x) 
is an approximate solution for the water-wave system under the following sense: 



Proposition 5.1. For any N gN, there exists 

N 

U''{t,x) = M{t,x) + V{t,x) = M{t,x) + J2^^^i^ ^^^h VieC'^{R+,H'' 

1=1 

and a small constant 6 = e~''°^ > such that for every I, one has the estimates 

mt)\E, < /iVCfc,,e-'^»(=^-^i)*, Vt > 

with some constant Ck^i- Moreover, U"" is an approximate solution of \1.(^) in the sense that 

dtU'' - T{U'') = Rap 

where Rap satisfies the estimate 

\Rap\E^ < C;v,./.'^<^^+^e-°(^+i)(--^^)*, for t > 0, 

In order to prove this proposition, we use again the Taylor expansion of J-" 

^ 1 
^(M + V) = F{M) + Y, j,D'HM]{V, ■■■,¥) + RnA^) 

1=1 

where the first derivative of J- is DJ^ = JA[M] where 

J _ f ^ l\ WAf] - /^~^A/ + 9 + ZMG[r]M]{ZM-) + ZmOxVm vmQx - ZmGItjm' 
" V-1 0) ' ''^'"J - V -dxivM-) - G[vm]{Zm-) G[mi] 
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with the notations 

Zm = Z[r]M,(pM], VM = v['nM,fM], Vm = 'P[riM,<fM]- 
The notations Z and v are introduced in Proposition 13.11 while V is defined by 

V[7],ip]u = 69, ((1 + {d^rjfr'ld^uy 
We shall also introduce the notations 

Zi = Z[Qi], Vi = vm, Vi=V[Qi], i = l, 2. 

Plugging V{t,x) = J2i=i^''^i{'^>^) i'^to the system leads to linear problems with source terms for 
Vk- The system for Vi is 

dtVi - JA[M]Vi = -Rm, (5.2) 

the system for V2 is 

dtV2 - JA[M]V2 = ^D^T[M]{VuVi) 
and the general equation for Vi is 

^ 1 

dtVi-jA[M]Vi = Y, E ziD'nMm„...,vi^). (5.3) 

p=2 i<li lp<l-i 

Before solving these systems, let us fix and recall some notations to be used in the remaining 
part of the paper. For U{t,x) = {Ui, U2Y, we define 

0<a+/3<fc ^* ^""^ 

\U{t)\y^, = sup \d^d^U{t,-)\L^, 

a+P<k 

where 

Note that X^ is the natural energy space for the water-waves system. 

5.1. The homogeneous linear system. The main ingredient in the proof of Proposition 15. II will 
be a rather precise estimate of the growth rate of the fundamental solution of the linear homogeneous 
equation 

dtV - Jk[M]V = (5.4) 

which corresponds to the linearization of the water-waves system about the multi-solitary wave M . 
As before, with 

1 0' 



^ ' -Zm 1 

we can perform the change of unknowns U = RV to get a simpler linear system which is equivalent 
to ([531) 

dtU -JL[M]U = {) (5.5) 

where 

-Vm + g + aM vmQx 



V -dx{vM-) G[riM\ 

is a self-adjoint operator with the notation ajv/ = a[M] = vmQxZm + QiZm and we will take 
Lm = L[M],Gm = G[r]M] for convenience. We shall also use the notations 

ai = a[Qi\, Li = L[Qi], i = 1, 2 
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and so on. Note that this is a short hand for: 

Liu{t,x) = L[Qc^{x — cit)]ii(t,x), L2u{t, x) = L[Qc2{x — h — C2t)]u{t,x). 

The main result of this section will be: 

Theorem 5.2. For e € (0,e*], there exists ho and C > such that for every h > Hq, the solution 
of (j5.5p with initial datum U{t) at t = t satisfies the estimate 

\U{t)\x>^+Y^ laf [72^2 < h-^Ck{\U{T)\x,+Y, |9rf/2(r)U2)(l+eo(c2-ci)(t-r)) V«(^^-^i)(*--)/2, 

a<k a<k 

Vt > r > 0. 

Remark 5.3. Note that in the above estimate, the right hand-side can be expressed in terms of 
usual Sobolev regularity of the initial data by using the system ()5.5p to express the time derivatives 
of the solution at the initial time t = t. Let us denote by SM{t,T) the fundamental solution of 
the (nan- autonomous) system (I5.5p . The estimate of Theorem 1 5. ^1 can thus be rewritten under the 
form: for every k < 0, there exists Ck{h, eo) > such that 

\SM{t,T)U\E, < Ckh-*{1 + eo(i - T)'=)|C/|^.we^»('^2-^i)(*-^)/2^ Vt > r > 0. (5.6) 

We do not need for our argument a sharp estimate of the number s{k) . A straightforward possibility, 
is to take s{k) = 2k + 1 since a time derivative of the solution always costs at most two space 
derivatives of the initial data. The meaning of this result is that when each solitary wave is stable 
(this corresponds to e small), then by choosing h sufficiently large, we can get an arbitrary slow 
exponential growth rate for the fundamental solution of (|5.5p . the special form of this growth rate 
that we have chosen is just one that it is sufficient for the proof of Proposition I5.il Note that 
the shape that we have chosen is linked in particular to the decay rate of the remainder Rm in 
Proposition \4-l\ 



We shall split the proof of the above estimate into many steps. For notational convenience, we 
shall give the proof only in the case r = which gives the worse constraint of /iq- The general case 
can be deduced from this one by replacing t by t — r, x by x — cit and thus in the multi-solitary 
wave M, h hy h = h + {c2 — ci)t > h. 

During the proof C is a positive number which change from line to line but which is independent 
oi h for h > 1 and t for t > 0. 

We shall first define a decomposition of unity in order to localize our energy estimates in the 
vicinity of each solitary wave. We take x^ ^ C°° (M) such that 



X'{x) 



1, x<0, 
0, X > 1 



and we define 



Xl{t,x) =X°( ^ —)^ Cr,^ = ^-^5-^, X2{t,x) = 1 -Xl{t,x). 
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Finally, we take 



Xiit,x) = ^ , , X2it,x) = ^ (5.7) 

{Xi + Xi)^ (XI + X2)2 



Note that these functions are smooth and bounded and defined such that xf + xi ~ ^■ 
The main properties of these functions that we shall use are: 
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Lemma 5.4. The above Xi satisfy 



y/3,\|3\>l,\^l,X^it,x)\<^, i = l,2. (5.8) 

Moreover, for any e > 0, we have 

|g-.|x-cit|^^| < ^^ |^_,|^_,,_e2i|^^| ^C^^ yt>0,x€R,h>l (5.9) 

for some C^ > 0. 

Proof. The estimate ()5.8p is clear. For the first one in ()5.9p . we observe that X2 is supported in 
X > Cmt + /i/4. Since in this region x — cit > {cm — ci)t + h/A > 0, we immediately get 

I -e|x-Cit| I < „-e(Cm-Cl)t -£-| < j: 

1"^ A^l r^ '^ "^ ^~^ h ' 

The second estimate follows by observing that xi is supported in x < h/2 + c^i- This completes 
the proof of Lemma 15.41 D 

5.2. Lower-order energy estimate. The first step in the proof of Theorem 15.21 will be to prove 
the estimate for k = 0. We shall thus prove 

Proposition 5.5. Under the assumptions of Theorem \5.SX we have the estimate: 

\U{t)\\, + \U2{t)\l, < c{h^U{T)\\, + |C/2(T)|i.)e^°(^^-^^)(*--)/^ Vt > r > 0. 

Proof. Again, we shall give the proof only for r = 0. 
Let us consider the energy functional 

Ei{U{t)) = {LmU, U) - cMxiU, xiU) - C2{Ax2U, X2U) 

with A the symmetric operator 

We shall prove the Ei is an almost conserved quantity with some positivity property thanks to 
Proposition 13.111 applied to each solitary wave. 
We first write 

^^^Ei{U{t)) = {dtU, LmU) + ^{[dt, Lm]U, U) - cMxidtU, xiU) 

- C2{Ax2dtU,X2U) - ci{A{dtXi)U, xiU) - C2{A{dtX2)U, X2U) 
:=/i+/2 + ---+/6. (5.10) 

We will estimate these terms one by one. Towards this, let us set 

C/i = xiU, and C/2 = X2U 
with ^,1^1 2 xf = 1 from (j5.7p . We shall use very often the following norm equivalence properties. 
Lemma 5.6. There exists C > such that for every h > 1, we have 

\U\j,o<C{Y,\Ufxo + lHD)U2\l2) 

i=l,2 

and also 



Y.\U'\xo<Cm'xo + l\<D)U2\l.] 



h' 

i=l,2 

with k{^) is a smooth cut-off function with k{£,) = 1 around ^^ = 0. 
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Proof. Let us prove the first estimate. From Lemma 15.41 we first easily get that 

i=l,2 "' 

For the estimate on U2, we use the basic commutator estimate 

\W.f]9\L^<\dJ\L^\g\L^. (5.11) 

Indeed, one can write [*P, f]g = (1 — d1)~^{dxfg) + [(1 — d1)~i,f]dxg- The L^ norm of the first 
term is clearly bounded by the right hand-side of (j5.1ip while the L? norm of the second can 
be bounded by C\dxf\L°=\g\L'2 by invoking [271 Proposition 3.6.B, estimate (3.6.35)]. This proves 
([5TT]1 . Estimate ([5TT]l yields 



|W2li2 < Y. \'^XIU2\1. < Y. (|[q3, X^]X^U2\l2 + |x.^(x.^2)li2) 
1=1,2 j=l,2 

< CYi\\X^U2\l. + mXiU2){] 



i=l,2 

The proof of the other estimate is similar. This ends the proof of Lemma 15.61 . D 

Now we can go back to the study of (I5.10p . First of all, from the linear system ()5.5p . we have 

h = {JLmU, LmU) = 0. 

Next, from the definition of Lm one can compute that 

ro ^ , _ f-[dt,'PM]+dtaM idtVM)dx\ 
^*' '''^"V -dAidtVM)-) [OuGm]) 
and so by combining this with the decomposition of unity leads to 

j=l,2 i=l,2 

where 
2/2/? = Y i([^*' ^Af - Li]x^U,XiU) + m, Lm], X^]XiU, U)} 

j=l,2 

= Y {((-[^t' 'Pm - Vi] + dt{aM - ai))xiUi, XiUi) + m{vM - Vi))^x{x^U2), X^Ul) 

i=l,2 

+i[dt,GM - Gi]xiU2, x^U2) - {[[du Vm], Xi]XiUu Ui) + m, GM],X^]XiU2, U2)}. 
For the remainder term I2R, we have by using Lemma 15.41 the estimate 

\l2R\<lc{\U{t)\j,o + \U2{t)\l,). 

Indeed, let us study for example the first term in more details. We have 

Vm = bdx 3- , 

\{l + \dxVM\T^J 

and thus 

r« T. T.1 i.F> ( l^if^ dxr]isdx{riM - rii) . \ 
[dt, Vm - Vi] = hdx\ / [5j, —ds]dx- 

\Jo {1 + \^xr]^s\^)^ J 

with the notation dxTjis = dxTji + sdx{riM — Vi)- Consequently, by using again Lemma E31 we get 

\[dt, Vm - Vi]x^Ul, XiUi)\ < lc\Ui\l^. 
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One can estimate the other terms in the definition of I2R by using the same arguments, in particular, 
for the terms involving the Dirichlet-Neumann operator, we also have 

Lemma 5.7. There exists a constant C such that the following commutator estimates hold 

{[dt, Gm - Gi]x^U2, X^U2) < -^^(1^2112 + \U2\I2) (5.12) 

m, Gm], X^]XiU2, U2) < lc{\W2\l2 + \U2\I2) (5.13) 

where Xi (^ = 1, 2) are defined in ^5.7\ ). 

This lemma will be proven in the appendix. 

For the I^ term in (|5.10p . we get by using the system (j5.5p that 
^3 = -ci{AxiJLmU, xiU) = ciidxixiLAiU), xiU) 

= -ci{LiXiU, d.,{xiU)) - ci{{Lm - L^)xiU, d^{xiU)) - ci([xi, Lm]U, d^{xiU)) 

with the remainder I^n defined by 

Im = -j^ci{[dx, Lm - Li]xiU, xiU) - ci([xi, Lm]U, dr,{xiU)) 

= -^'^lii^x, Vm -Vi- OM + ailxif^i, Xif^i) + ci((<9x^'A/ - dxVi)xiUi, dx{xiU2)) 

+ \ci{[d^, Gm - Gi]xiU2, X1U2) - ci([xi, Lm]U, d^{xiU)). 

Note that the structure of /3/j is very similar to the one of the I2R term above (basically dt is 
replaced by dx in the commutators) and hence by using the same arguments as above, we get 

\hR\<lc{mt)\'^o + \u2{t)\i2). 

In a symmetric way, we also have for the /4 term in (j5.10p that 

Im = ^C2{[dx, L2]x2U, X2U) + hn 



with 



\hR\<lc{\U{t)\j,o + \U2{t)\l2] 



Since the solitary waves have the dependence 

Ql = Qlix - Cit), Q2 = Q2ix - h- C2t), 

we have that 

[dt, Li] = -Ci[dx, Li]. 
This yields the crucial cancellation 

h + I3 + h = hR + Hr + Iar 
and hence, we obtain the estimate 

\h + h + h\<\c{\u{t)\l, + \U2{t)\l2). 

By using integration by parts and (|5.8p . we also easily get that 

\h\ + \h\<^{\Ui\H^+\U2\L^). 
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Summing up the estimates, we get from (j5.10p that 

\j^E^{U{t)) < \c{\U(t)\\, + \U2{t)\l.). (5.14) 

The next step will be to get a minoration of Ei{U{t)). By using the decomposition of unity (j5.7p 
again, one has 

E^{U{t)) = {LmU, U) - ciiAxiU, xiU) - c2{Ax2U, X2U) 

= Yl (^MXiU, XiU) - ciiAxiU, xiU) - C2{Ax2U, X2U) + Y, d^M, Xi]XiU, U) 

i=l,2 i=l,2 

:= Ih+Ih + Ih (5.15) 

where 

111 = {LmXiU, XiU) - ci{AxiU, xiU), 

112 = {LMX2U, X2U) - C2{Ax2U, X2U) and Ih = ^ ^^^M, X^]XiU, U). 

4=1,2 

We shall first handle Hi. We note that 

Ih = {LiXiU, xiU) - ciiAxiU, xiU) + {{Lm - Li)xiU, xiU) 

= (LiXiU, xiU) - {{Lm - Li)xiU, xiU) (5.16) 

where 

-Vi + g + vidxZi + dtZi {vi - ci)da: 



^'"' -dA{vi-ciy) G[m] 

Noticing that Qi = Qi{x — cit) and so dtZi = —cidxZi, we can rewrite Li as 

^ _ f-Vi + g + {vi- ci)dxZi {vi - ci)dx 



-d.{{vi-ci)-) G[7?i] 

Note that the operators Li is the same operator as Lc^ studied in section [3^ except that it coeffi- 
cients depends on Qi = Qi{x — cit), we have 

TcitLi = LciTcit 

where Txq is the translation operator 

{%oU){x) = U{x + xo). 

Since Tct is an isometry on L^ and X^, Proposition 13.61 applies to Li. Let us use the notation that 
Q'i{x) = dxQci{x) and define 

UHt,y) = U\t,y + cit) = {Tc,tU^){t,y) = xi{t,y + cit)U{t,y + cit). 

Thanks to Proposition 13.11] we can use the decomposition 

U\t,y) = ai{t)JRiQ[{y) + f3i{t)RiQ[{y) + W\t,y) (5.17) 

such that W^ satisfies 

{W\{i]'i,oy)=0, {W\jRiQ[) = 0. 
Note that we use again the short-hand Ri = R^. Since by definitions, we have 

{LiU\U^) = {T-c,tLc,Tc,tU\U^) = {L,,Tc,tU\Tc,tU^) = {L,,U\U^), 
we get from Proposition 13. Ill that 

(Li^S [/^)>co|t^^|^o-C|ai|2. 
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Moreover, as in the estimate for I2R above, we also get from Lemma 15.41 that 

|((Lm - Li)xiU, xiU)\ < ^Crnxo + \U\h)- 
Consequently, we get from the two previous estimates and (j5.16p that 

Ih > co\W'\j,o - C|ai|2 - ^C{\U\lo + \U\l,). 
In a symmetric way, we can set 

U\t, y) = {Tc,t+hU^){t, y) = U^{t, y + h + cgt) = X2(t, y + h + C2t)U{t, y + h + cai) 
and use the decomposition 

U\t,y) = a2{t)JR2Q'2{y) + Ht)R2Q'2{y) + W\t,y) (5.18) 

with 

(^', (r?2,0)*)=0, {W\JR2Q'2)=^ 

to obtain that 

Ih > co\W^\^^o - C|Q2p - lc{\U\j,o + \U\l,). 
Moreover, we also have from Lemma 15.41 the estimate 

\n3\<lc{\u\j,o + \U2\l.) 

by using again that the commutator always involves at least one derivative of Xi- 
In view of the decomposition ()5.15p , we have thus obtained that 

Ei{U{t)) > co{\W%o + \W^\xo) - C{\ai\^ + |a2p) - lc{\U\j,o + \U2\I2). (5.19) 

In order to conclude, we still need to estimate |C/2|l2, |aj| and \f3i\{i = 1,2). 

For the L^ norm, let us choose k{D) where k G Cq°(M) and k(^) = 1 around ^ = 0. From the 
linear system ()5.5p . we have 



dtU2 = (Pm - cimWi - gUi - VMd^U2, 
and thus we obtain 

^j^\k{D)U2\12 = {K{D)dtU2, k{D)U2) 

= {k{D){Pm - aM)Ui, k{D)U2) - g{K{D)Ui, k{D)U2) 

-{K{D){vMd^U2), k{D)U2). 

By using that k is compactly supported, this yields 

]^^^\K{D)U2\l2 < C{\Ui\h. + |W2|lO {\W2\L- + \<D)U2\r2) 
and hence, we obtain from the Young inequality that 

^^^\K{D)U2\l2 < C(6- Vixo + e\<D)U2\l2) (5.20) 

where e is a small constant to be fixed later. 
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To estimate Oj, we use the decompositions (|5.17p . (jS.lSp of C/* {i = 1,2). We have 
|«i = \jRQ^\2^ 9tixim{t),JR,Q[) 



1 



, ^ {{dtXi{t))U{t), JRiQ'i) + iximJLMU + cid.,U){t), JRiQ[) 



\Tn\n2 [{{dtXi{t))U{t), JRiQ[) + {xi{t)J{L^ -L^){t), JRiQ[] 



+{xi{t){JLiU + cid^U){t), JRiQ[] 
with the notation that f{t,x) = f{t,x + cit). We also have 



{xi{t){JLJJ + cAum, JRiQ'i) 

= {xi{t)J{L^-ciJdJJ){t), JRiQ[) = ixiit)JLi[Qi{y)]Uit), JR^Q\) 

= ([xi(t), JU[Qi{y)]]U{t), JR,Q[) - iximit), U[Qi(:y)]J^RiQ'^) 

= ([xi(t), jZi[Qi(y)]]C7(t), JiJiQl). 
Indeed, to pass from the second to the third hne, we have used the crucial cancellation 

(Li - ciJd,)J^RiQ[ = -LiRiQ[ = -L,,R,,Q'^^ = 0. 
By using again Lemma 15.41 to estimate the other terms, we thus obtain 

\j^ai\<lc{\U\xo + \U2\L^). (5.21) 

In a symmetric way, we also get for \a2\ that 

\f^a2\<lc{\U\xo + \U2\L^). (5.22) 

We still need the estimates of |/3i| and |/32|. By using ()3.40p . we note that we can write 

Pi=Pi-ai{JRiQ\,{7][,0Y) 
where 



/3i 



In particular, thanks to (j5.2ip . we obtain that 



lr?'|2 



||/3i| < ||a| + lc{\U\xo + \U2\l^). (5.23) 



To estimate /?i, we directly compute 
d ~ , . 1 



i/3i(t) = j^dtiUHt), {rj[,OY] 



hi I 



l^({{dtxi{mit), {v[,0Y) + {ximdtU + cid,U){t), (7/1,0)*; 



L2 



1 



lL2 



I ,|2 v((5iXi(i))f^(t), (Vi,0)*) + (xi(t)Ji:Mi7(i), (r?l,0)*^ 



+ci(xi(t)a.i7(t), (r?;,0)*) 
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where in particular 



{xi{t)JLMUit), (77;, 0)*) = {xi{t)J{LMU - LiUrn, (r/'i,0)*) + ([xi(t), JLi]U, {^[,0]') 
and 



ci(xi(tR^(t), {v[,oy) = -ci{{d,xim{t), {v[,oy) + ci{d,m{t), {v[,oy). 

Therefore, we get 

|^/3i(t)| < C\{JL,Tj\{r^[,0))\ + ^C{\U\xo + \U2\L^). 

To conclude, we use the decomposition (I5.17p . since LiR\Q'^ = 0, we have 

|jZiI7\(??;,0))| <C{\ai\ + \W^\x<^) 
and hence, we get from (I5.23|) that 



|^/3i(t)| < c{^{\U\xo + \U2\l^) + |ai| + iW^^lxo). (5.24) 



Similarly we have 

\j/2{t)\ < C {^{\U\x^ + |C/2|lO + l«2| + iW^'lxo) . (5.25) 

From dLZU, ([522]), ([521]), dS^S]), we finally obtain 

^^ai?<C\ai\\{\U\x^^ + \U2\L^), i = l,2 (5.26) 

at h 

1 1 

^|A|'<C|A|(-(|[/|xo + |t/2|L2) + |a,| + |t^^|xo), 2 = 1,2. (5.27) 

To combine (I5.26P , (I5.27P , (I5.14p , (I5.20p and (I5.19P in an efficient way, we define a weighted energy 
Ei{U{t)) = lh-2E,{U{t)) + ^(|/3i(t)P + |/32(t)P) + ChH\a,{t)\' + \a2{t)\') (5.28) 

From ([STip and (fOOD . (fOGl) . (fOTD . we obtain 

|^i(^(t)) 

< /i-3C[/i-i(|[/||o + \K{D)U2\l2) + /i-3(|/3i| + \p2\){\U\xo + \U2\L2) 

+hH\l3i\ + l/32|)(|ai| + lasl + |T^^|xo + \W^\xo) + /i^3(|ai| + \a2\){\U\xo + IC/sIlO] 

< ^-3C[/i5(|Vl^i|2^o + \W^\lo) + (|/3iP + |/32p) + /i^(|ai|2 + |a2p) + /i-^lC/glia)] . 

To get the last line, we have used the decompositions of U^, f/^, (|5.17p . (j5.18p and Lemma [^ 
Now, let us define 

F{t) = hk\W'\lo + \W^\Jco) + (l/SiP + l/32p) + h'mail^ + |a2p) 
and integrate the estimate for Ei{U{t)) with respect to time from to t, we obtain 

Ei{U{t)) < Ei{U{0)) + h~^C / (F(s) + h~^\K{D)U2{s)\l2)ds, 

Jo 
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by using that |f/2|L2 ^ C{\i^{D)U2\l2 + \U\x«)- On the other hand, from the definition of Ei(U) in 
(j5.28p and the estimate (j5.19p . we also have 

EiiUit)) > lhho{\W%o + \W^\j,o) + lm\^ + m^) 

+hho{\ai\^ + |a2p) - h~-2Ci ^ |C/*||o + \k{D)U2\12) 

1=1,2 

> coF{t)-h^ic\K{D)U2\l2 
for some cq > 0. This yields 

F{t)-h-'^C\K{D)U2{t)\l2 <h^C\Umlo+h-'^C / {F{s) + h-'^\K{D)U2{s)\l2)ds. (5.29) 

Jo 

By taking the integral of ()5.20p with respect to time, we also have 

\K{D)U2it)\l2{t) < \U2mh +C I ie-^F{s) + €\fiiD)U2is)\l2)ds (5.30) 

Jo 

where the constant e will be fixed later. We take the sum of (j5.29p and (j5.30p multiplied with some 
A > that will be also chosen later to get that 

F{t) + {X-h-"2C)\KiD)U2it)\l2 



1 '■^ 



<h^C\U{0)\j^o + X\U2mi2+Ch~-^ / {F{s) + h—i\K{D)U2{s)\i2)ds 

Jo 

+XC I {e^^F{s) + e\K{D)U2{s)\l2)ds. 
Jo 

Let us choose A and e in order to satisfy the following conditions 

A - h-^C > A/2, C/i-3 < 7/4, /j-i < A/2, ACe~^ < 7/4 and e^ < A/2, 

where 7 = eo(c2 — ci) comes from Proposition 14.11 One can choose for example A = g^^, ^ = 2 
and h large enough to get our final energy estimate 

F{t) + ^\n{D)U2{t)\i2 < h'^cmorxo + x\u2mh + 1 l\ns) + ^\K{D)U2{s)\i2)ds. 

Applying Gronwall's inequality, we get 

F{t) + ^\K{D)U2{t)\l2 < (/i^C|C/(0)||o + X\U2mh)e^'^\ for any t > 0. 

From the definition of F{t) and Lemma |5.6|, we get that there exist constants cq and cq such that 

F{t) = hk\W'\%, + \W^\j,o) + m\^ + m^) + hk\ai\^ + \a2\') 

By choosing h such that ^ < ^, we find 



> CO 5^ \W{t)\lo > (1 - ^)c-o|C/(t)|^o - "f\^{D)U2{t)\l2. 

i=l,2 



F(t)>ico|C/||o-^|«:(i^)C/2li2. 
Summing up the estimates above, we finally obtain 

y It^lxo + ^I^(^)f^2li2 < {h'2C\U{0)\j,o + A|C/2(0)|i.)e^''('=2-'^l)*/^ for any t > 0, 
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that is to say 

|t/|^o + \U2\I2 < C{hl\Umxo + |f/2(0)|i2)e^o(^^-^l)*/^ t > 0. 

This ends the proof of Proposition 15.51 

5.3. Proof of Theorem 15.21 We shah prove Theorem 15.21 by induction. The k = case was 
aheady obtained in Proposition [53J For the sake of clarity, before the general induction argument, 
we shall explain the proof for the k = 1 energy estimate. 

Order- 1 energy estimates. Note that since the coefficients in the linear system ()5.5p depend on 
t and X, neither dt nor dx has a nice commutation property with the equation. We shall thus use 
the operator 

D{d) = dt + cixldx + C2xldx (5.31) 

in order to take derivatives of the equation. This will take into account the fact that the solitary 
waves depend on x — Cjt. 

By applying this operator on both sides of system (15. 5p . we obtain for D{d)U the system: 

dtD{d)U - JLMD{d)U = [du D{d)]U - J[Lm, D{d)]U. (5.32) 

A further computation shows that 

[dt, D{d)]U = ci{dtxl)dxU + C2{dtxl)dxU, 

and 

J[Lm, D{d)]U =J[Lm, dtp + JY1 c*[^m, xldx]U 

i=l,2 

=J Y, [xULm - Li, dtp + x}[U, dtp + Ci[LM, xhdxU + c^x^Lm - U, d^p 

1=1,2 

+ x![Li, c,dxp) 

=J Y, {x1[Lm - Li, dtp + c,[Lm, xl]dxU + c^x^Lm - L,, d^p) := JSU. 

i=l,2 

(5.33) 
The crucial fact that we have used in the above computation is the cancellation: 

J J2 XiiLi, dtp +jY1 x-[Li, cAp = jY. ^^[^^' ^* + ^*^-]^ = 

i=l,2 i=l,2 2=1,2 

We can thus rewrite the system (j5.32p for D{dp as 

dtD{dp = JLMD{dp - JSU + [dt, D{d)p := JLMD{dp + Fip) (5.34) 

where we will take Fip) = —JSU + [dt, D{d)p as the source term. 

From Proposition 15.51 we have for the fundamental solution SM{t,T) of the linear equation (|5.5p 
the estimate 

\\SM{t,T)[[xonL^^xonL^ < /l3Ce^o(^^-^l)(*-)/^ Vt > r > 0. 

Consequently, by using Duhamel's Formula, we can rewrite ()5.32p as 

D{dp{t) = SMit,0){D{dp){0) + [ SM{t,T)F^{T)dT, 

Jo 
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so we get 

\D{d)U{t)\x^ + \D{d)U2\L^ 

a=0,l -^0 

(5.35) 

We still need to estimate the source term Fi . We shall first use the elliptic regularity of the leading 
spatial operators in (j5.5p to prove that time derivatives control higher order space derivatives. 
For any norm || • || on x dependent vectors, we use the notation 

\mfu\\= Y. \\d\u\\. 

0<l<k 

Lemma 5.8. Any smooth solution of i5.5\) satisfies the following a priori estimates: 

V/ > 0, m > 0, 3C,,i, \di{U,,U2)\^^^.^^^^, < Q,J{dty-^HUi,U2)\^^^,^^^^i. (5.36) 



Proof. This is an elliptic regularity result. Let us start with the case / = 0. We first rewrite ()5.5p 
under the form: 

Gm U2 = dtUi + dxivMUi), , . 

VmUi = dtU2 + {au + g)Ui + VMd.U2. ^^-""'^ 

The operator Vm is an elliptic operator of order two, therefore, by classical elliptic regularity results, 
we immediately obtain from the second line of the above system 

and thus by the interpolation inequality 

we actually obtain 

|C/i|^™+5 <C„.{\dtU2\^^^i+\U2L^ + \Ui\L2). (5.38) 

In a similar way, since the surface rjM is smooth, the Dirichlet-Neumann operator G[?7jv/] is an 
elliptic operator of order one. Actually, we have 

G[7]M] = \Da:\+R{t,X,D^,) (5.39) 

where i? is a pseudo differential operator of order zero. We refer to [T7], [28] for the proof. Note 
that we are in dimension one thus the principal symbol is very simple. Consequently, we also get 
from the first equation of (j5.37p 



\U2\h^+'^ < Cm(\dtUi\Hn.+l + \Ui\Hrn+2 + |^2|l2J- (5.40) 

By using again that 

l^lb.n+2 < 6\Ui\^„^^S+Cs\Ui\l2, IC/2L+3 < 6\U2\h^+2+Cs\U2\l2, 

we get from (|5.38p . (j5.40p by choosing 6 sufficiently small that 



|C/l|^„,+ 5 + \U2\h^+^ < C„^[\^tUl\Hm+l + \dtU2\^^+i + \U\l2J. 

We have thus proven (j5.36p for / = 0. 

We then proceed by induction on /. Assume that the result is proven for / — 1 time derivatives. 
Let us apply 9| to (|5.37p . From the second line, we get 

VudiUi = d[+^U2 + d[{{aM + 9)Ui) + d[{vMd,U2) - [dlVMpi := F^ 
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and we observe that the right hand-side satisfies the estimate 

iF\„+i < Q,m{\di+'u2\^^^i + i(5t)'c/2i^„+3 + my-'uii^^,^, 

and thus we get by ehiptic regularity that 

\diUi\^„,^s < Q,„(|a^iC/2|^„+i + K5t)'C/2|^„+3 + \{dty-'Ui\^,„^s). (5.41) 

In a similar way, for the first line of (j5.37p . we get 

GMdiU2 = dl+^Ui + did^A^MUi) - [di, Gm]U2 ■.= f\ 

By using Proposition 13.11 (6) to compute the commutator [dl, Gj\/][/2j we obtain for the right hand 
side the estimate 

Consequently, from the ellipticity of the Dirichlet-Neumann operator, we also get 

\diU2\Hm+2 < CLm{\dl^^U^\H^+i + \{dtyUi\Hn.+2 + \{dty~^U2\Hm+2y (5.42) 

By combining (|5.41|) and (j5.42|) we get 

\dliUi,U2)\^„^,s^^^^, < Q,^(|5^i(C/i,C/2)|^„,,^^.,i 

+ |9i(^i,C/2)|^„+,^^„.^3+|(9,/-HC/i,C/2)|^„^5^^„^,). 
To conclude, it suffices to use the interpolation inequality 

|9i(c/i,c/2)|^„,,^^.,3 < 5|ai(c/i,c/2)|^^,5^^„^, + c,|ai^U. 

and the induction assumption. This ends the proof of Lemma 15.81 D 

As a consequence of Lemma 15.81 we get the following inequalities that we will use many times: 

• By using the lemma with / = 0, ??t, = 0, we have 

19x^1 3 , < C(\{dt)U\ , i) 
and hence, since, we can use again that 

\d^U\ , i<e\d^U\ 3 +CAU\l2 (5.43) 

we get that 

|5.f/Uo < e(|5i[/Uo + \dtU2\L^) + C,{\U\xo + \U\l2) (5.44) 

for any e > 0. 

• We can also use Lemma 15.81 with / = q, 77i = /3 — 1 for any a and /3 such that a + (3 = k 
and /3 > 1. We obtain 

and we can iterate the process to obtain 

\d?d^^U\ 3 <\(dt)''U\ g.3 ,^< Ck\{dt)^U\ , i. 
Thanks to ()5.43p . we thus obtain 

laf 9f C/^o < e{\dkU\xo + \d^U2\L^) + a{\U\x.^. + \{dt)'~'U2\L^) (5.45) 

for a + /3 = /t, /3 > 1. 
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We can come back to the estimate for the source term Fi{U{t)) in the right hand side of (j5.35p . 
By using the expression of Fi{U{t)) in (|5.33p . we get 

\FiiU)\xo < \JSU\xo + \[dt,Did)]U\xo 

< \J Yl ^fi^M - Li,dt]U\xo + \JY1 <^iiLM, xf]dxU\xo 

i=l,2 4=1,2 

+1^ Yl ^^xI[Lm - U, d^]u\xo + \[dt, D{d)]u\xo. 

i=l,2 

By using again Lemma 15.71 and Lemma 15.41 as before, we obtain 

\FiiU)\xo < ic(|(5,.[/2, d'',UiY\xo + \dxU\xo + \U\xo), 
and hence, by using Lemma l5.8( we find 

\Fi{U)\xO < ]-C{\dtU\xO + \U\xo + \dtU2\L^ + \U2\l^). 

h 



In a similar way, we also have 

\Fi{U)\l2 < l-C{\dtU\xo + \U\xo + \dtU2\L^ + \U2\l2). 
h 

Going back to (j5.35p . we obtain 

\D{d)U{t)\xo + \D{d)U2{t)\L2 

Jo 
It remains to deduce an estimate of the time and space derivatives from the previous estimate that 
gives only a control of the D(d) derivatives of U. We first observe that 

\Did)Uit)\xo + \Did)U2{t)\L2 

> \dtU{t)\xo + \dtU2{t)\L2 - Y Ci{\xldxU{t)\xo + \xfd.U2{t)\L2) 

i=l,2 

> \dtU{t)\xo + \dtU2{t)\L2 - C{\dxU\xo + \d,U\L2). 
We know from (15.44P that 

\Did)Uit)\xo + \Did)U2{t)\L2 > \dtUit)\xO + \dtU2{t)\L2 - Ci\Uit)\xo + \U2{t)\L2). 

We have thus proven that 

\dtUit)\xo + \dtU2{t)\L2 - C{\U{t)\xo + \U2{t)\L2) 

< /iiCe^o(c.-cOt/4(|^(o)|^^ + 1^^(0)1^^ + \dtU2mL2) (5^4g) 

+ /i-lc [ e^«(^^-^i)(*-^)/^(|ai[/(r)Uo + |C/(r)Uo + \U2{t)\l2 + \dtU2{T)\L2)dT. 
Jo 

From Proposition 15.51 we have 

\UiT)\xo + \U2{t)\l2 < /i2C(|^(0)Uo + |C/2(0)U2)e^»('=^-'=lW^ for r G [0, t], 

and therefore 

h-ic [ e^o(^2"^i)(*"^)/4(|C/(T)|xo + \U2iT)\L2)dT <th--2C{\UiO)\xo + |f/2(0)|i2)e^o(^^-^^)*/^ 
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consequently, we finally deduce from (|5.46p that 



\dtU{t)\xo + \dtU2{t)\L2 < h-^C{l + h-h)e'^^''-''Wi{\U{0)\xi + |C/2(0)|i2 + \dtU2mL2] 







Note that we could avoid the additional algebraic growth but that we do not need to refine. 
The Gronwall's inequality yields that there exists a constant Ci such that 

\dtU{t)\xo + \dtU2{t)\L2 < /i^C7i(l + /i-ii)e^o('=^-'=i)*/^+'^"*^i*(|?7(0)Ui + Yl 15^^2(0)1^2). 

a=0,l 

This yields the desired estimate for the time derivative by taking h large enough. We finally deduce 
from (15.44P that the same estimate also holds for \dxU{t)\xo- This ends the proof of the order-1 
energy estimate. 

Higher-order energy estimates. We will do this by an induction argument. First of all, assume 
that we already have the estimates for |C/|jffe-i and |K(D)9f [721^2 with k > 2 and a < k — 1: 

a<k~l 

a<k-l 

As previously, we get start with the estimate of \D {d)U\xo- One can write the system solved by 
D''{d)U as 

dtD\d)U = JLMD{dfU + Fk{U) (5.47) 

where 

fc-i 
Fk{U) = YD{dy{J[D{d), Lm] + [dt, D{d)])D{df~^-'U 

i=0 

will be considered as the source term. 
From the Duhamel formula, we find 

D{d)''U{t) = SM{t,0)D{d)''U{0)+ [ SM{t,T)FkiU{T))dT 

Jo 

and therefore, we again obtain from Proposition 15.51 that 

\D{d)'Uit)\xonL^ 
< /i3Ce^o('=2-'=i)*/4 (1^7(0) l^fc + V |9fC/2(0)|i2) + /i3C /" e^"(^2-^i)(*-^)/^|Ffc(r)|xonL2(ir. 

a<k-l -^0 

(5.48) 
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It remains to estimate the source term Fk{t). We have 

\FkiU)\xo 



k-l 

< Y^ {\JD{dy[D{d), LM]D{d)'^''-^U\xo + \D{dy[dt, D{d)]D{d)'^-'-^U\xo) 

1=0 

k-l 

< E E [\JD{dyx][LM - Lj,dt]D{d)'^-^-'U\xo + c,\JD{dy[LM, x]]d.D{d)''~^-'U\xo 

i=0 j=l,2 

k-l 

+c,\JD{dyx][LM-L„ d,]D{d)''-'-'U\xo]+Y,\D{dy[dt, D{dyD{d)''-'-'U\xo. 

i=0 

Consequently, by using again the same arguments to estimate the commutators and the observation 
(already used in the order-1 estimate) that the commutator [D{d), Lm] = S (computed in (|5.33|) ). 
acts like j^Lm, we can get that 

\Fk{U)\xo + \F,{U)\l2 < \Ck E (l(9f Sf ^2, af 9f+iC/i)*|^o + \dfdl,U\xo) + \c\U\x.-i, 

a + /9 = fc 

a<k-l 

and so by using Lemma 15.81 again we arrive at 

\Fk{U)\xo + \Fi,{U)\l2 < lCk{\d^U\xo + \U\x^-^ + E \drU2\L^). 

a<k 

Going back to (j5.48p one has 

\D{d)''U{tyxo + \D{d)'^U2{t)\L2 

< /iiCe^o(-^-^)*/4(|C/(0)|^. + Y, \drU2m\L-)) 

a<k 

+h-iCk re^o(^^-^i)(*--)/M5f^(T)|^o + |C/(r)|^._, + ^ |5r^2(r)|L2)cir. 

•^0 a<k 

For the left-hand side, we observe that 

\D{d)''U{t)\xo + \D{d)''U2{tyL^ > \d^U{tyxo + \d^U2{tyL--C{ Y \drd!u{t)\xo + \U{t)\x.-i), 

a+/3<fe 

a<k-l 

which together with (J5.45P allows to get 

\D{d)^U{t)\xo + \Did)'^U2{t)\L2 > \d'^U{t)\xo + |9f C/2|i2 - C7(|[/(t)|^.-i + ^ l^f ^^sIl^). 

a<k-l 

Consequently, by using the induction assumption to estimate |C/|jffe-i + J2a<k-i l^f ^2|l2, we find 

\d'^U{t)\xo + \d^U2\L^ < /i^Cfc_i(l + /i-!t)'=-ie^o('=^-^^)*/4+^"^^i*(|[/(o)|^. + Y \d?U2m\L^) 

a<k 

+h-lCk I e^«(^2-'=i)(*-^)/4(|9f[/(r)|^o + |a,^'f/2(r)|i2)dr. 

JO 

Using the Gronwall's inequality as before yields 

\dtU{tyx^ + \d^U2\L2 < /i3Cfe(|C/(0)|^fe + 5^5f |[/2(0)|)(1 + /i-!t)Vo('=2-=i)«/4+/.-4cit_ 

a<k 
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^_^^ 3 

Finally, by using again Lemma 15.81 and by taking h large enough such that h^^Ci < eo(c2 — ci)/4 
we conclude that 

\Uit)\x^ + E \9?U2\l- < h^CkiWm^u + Y, |5rf^2(0)|)(l + 6o(c2 - ci)tfe^'^-^--'^"\ Vt > 0. 

a<.k a<k 

This completes the proof of Theorem 15.21 

Remark 5.9. Note that in the proof of Theorem \5.^ the only information about the positions of 
the solitary waves are in the interaction estimates given in Lemma \5.4\ These localization estimates 
are true on any interval of time for which the distance between the centers of the solitary waves is 
bigger than h which is considered as a large parameter. We can use this remark to get an estimate 
for SMit,T) for < t < T. Indeed by using again the reversibility symmetry of the water waves 
system and the symmetries of the solitary waves of Theorem \1.1[ consider U{t,x) the solution of 
(j5.5p with initial data at t = t, then 

U{t,x) = {Ui{t - t, -x), -U2{t - t, -x)y 
is still a solution on [0, r] of (j5.5p with M{t,x) = Qi{x — cit) + Q2{x — C2t — h) replaced by 

M{t, X) = Ql{x + CiT - Cit) + Q2{X + C2T + h - C2t) 

and with initial data for U at t = 0. Consequently, we see that on [0, r], the center of the solitary 
waves are located at x = xi = —ci\t — t\ and x = X2 = —C2\t — t\ — h. Consequently, the slow 
solitary wave Qi is now located on the right. Nevertheless, we observe that for t G [0,t], we still 
have 

a^i — 3^2 > (C2 — Cl)\t — t\ + h > h 
and thus the solitary waves are still at least at distance h uniformly in t. Consequently, we still get 
as in Theorem \5.S\ that 

\U{t)\E^- </i^Cfc|C7(0)|^.w(l + eo(c2-ci)t'=)e'^, Vt G [0,t]. 
This yields in the original time variable 

\U{t)\E, < h\Ck\U{T)\Hsi,) (1 + eo(c2 - ci)(t - t)'=)e'^^, Vt G [0, r]. 

By combining with Theorem \5. 21 we thus obtain that the fundamental solution 5M(t, r) of system 
(j5.5p enjoys the estimate 

\SM{t,T)U\E^ </i^Cfe|C/|j^.w(l + eo(c2-ci)|t-r|*^))e'^^, Vt, r > 0. (5.49) 

5.4. Proof of Proposition I5.lt Construction of the approximate solution. Now we can 

go back to the study of the linear systems solved by V; (1 < Z < N) (j5.2p . ()5.3p in order to prove 
Proposition 15. 1[ 

We first note from the fact that (j5.4p and (j5.5p are equivalent via the transformation U{t) = 
RV{t), with R invertible, we get from Theorem 15.21 and (j5.6p 

|5j^(i,r)y|^. </iiCfc(eo)|F|H.w(l + |t-rhe^«(^2-^i)l*-^l/2, Vt, r > (5.50) 

where S^j{t,T) is the fundamental solution of the system (j5.4p . 

Let us go back to the construction of the approximate solution V{t,x) = X^;=i (5'V;(t,x) with 
S = e~^°^. For Vi(t), we have to solve 

dtVi - JA[M]Vi = -Rm, 

where the right hand side satisfies (see 15. ip the estimate 

\Rm\e^ < Cke-'"^'^-"'^', Vt > 0, A; G N. 
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We choose the solution 



Stiit,r)RM{T)dT. 
From the estimate (j5.50p of the fundamental solution, Vi is well-defined and satisfies the estimate 

1=0 

< /i3Cfc_i(eo)e-^o(^2-ci)t^ fQ^ ^>Q^ 
For the general case of T^ (2 < / < A^), we have to solve 

' 1 

dtVi-JK[M]Vi = Y, Y. -D'PT[M]{Vi,,...,Vi^):=Ri{t,x). (5.51) 

p=2 i<ij^,...,ip<]V 

We shall use an induction argument. Let us assume that we already have already built Vj, 1 < J < 
I — 1 that satisfy the estimate 

\VMe^ < /i^aj(eo)e-^'^°(^^-^^)*, Vt > 0. 
We get for the right-hand side of (j5.5ip that 

' 1 

\Ri{t,x)\E. = IE E ^DPT[M]iVi„...,Vi^)\ 

l<ll,...,lp<N 
ll+-+lp=l 

Taking 

/oo 
S^j{t,T)Ri{T)dT 

as a solution, we get thanks to (|5.5U|) 

2/-1 r°° 



2'-l ~ ' ^ -^leo{c2-ci)t 



This ends the proof of Proposition 15. li D 

6. The nonlinear problem 

After the study of the approximate solution f/" of water-wave system, we need to consider the 
remainder solution U = U — U°' where U is the solution of the water-wave system (II. 6p and U 

C Q -J- -1 GTI OC 

dtU^ = F{U'^ + U^) - F{U'') - Rap, t > 0, , . 

U^{0) to be fixed later '^ ' ' 

Proposition 6.1. Let m > 2, U°- £ W^^. and Rap G ^\l^'^)- There exists a solution U^ = 
{ri^, ip^Y e L°°([0,oo), H'^+^ X F'^+i) for ^UJH with a fixed initial value U^{0) such that 

H - Wv^Wl-- - I^^IIl- >o 

44 



and 

\U''{t)\^^^,^^^^r < C;v,^/i^^<^^+ie-(^+i)^«(^-^^)*, for any t G [0,oo). 

Proof. The proof is left to the reader, it suffices to use the same arguments as in the proof of 
Theorem II. 3[ D 

End of the proof of Theorem \1.SX Since we have aheady shown the global existence of the remain- 
der solution U^ , we know that there exits a (semi-) global solution f7(t, x) = U"'{t, x) + U^(t, x) = 
M{t, x) + V{t, x) + U^{t, x) to the water-wave system (jl.Sp . It only remains to describe the asymp- 
totic behavior of U{t) when t tends to -|-oo. Since 

N 

U{t) = M{t, x) + Y, ^^Viit) + U^it) 
1=1 

with 6 = e~'''' and with Vi and U that satisfy the estimates from Proposition 15.11 and Proposi- 
tion 16.11 we have 

This gives in particular 

\Viit)\Hs < h'-^CN^se-'^^^''-"'^', iW'mHs < /i'^C^,.5^+ie-(^+i)^o(^^-^i)* 
for any t G [0, oo). We thus get that 

lim \U{t)-M{t)\H- =0. 

i— >-+oo 

This ends the proof of Theorem 11.21 D 

7. Appendix: Proof of Lemma [5771 

7.1. Proof of (|5.1'2|) . We shall only prove the case i = 1, the case i = 2 can then be obtained by 
symmetric arguments. Let us recall that as in (13. 2p . (jS.ip . the Dirichlet-Neumann operator G[r]] 
can be defined by: 

G[ri]u = d^u%=o 

where u satisfies the elliptic system on the flat strip 5 = M x [—1,0] 

' yx,z-PVx,zu'' = 0, in cS 

with the notations 



and 9„ = n • P'Vx,z where n = — e^ is the outward unit normal to the boundary z = —1. Note that 
aP\ i_k 

^n \z=-l — H+T)^^- 

One can see by the Green's Formula that 



{G\n\u,v) = I pv,,,u^ • v,,y 



and so we have 



{[dt, G[rj\]u, u)= f dtPVx,zu' ■ Vx,,u' - 2 / PV^Aidtuf - Otu') ■ V^^.u' 
JS Js 

In the following we shall use the notations (j = M, 1, 2) 

G[vj] = G„ P[r],] = P, and u%]=u';. 
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We also set u = X1U2 for the sake of convenience. Then we can write 
{[dt, Gm - Gi]u, u) 

= I dt{PM - Pi)'^x,zu\i ■ ^x,zuii + / dtPiVccAuii - 4) ■ ^x,zu\^ 
Js Js 

+ I dtPiVx,zu\ ■ VxA^i - "?) - 2 /" (Pm - Pi)yx,zu\i ■ VxAidtu)ij - dtu\,j) 

-2 j PiV:,Aidtu)\i - dtu\i - {dtu)\ + dtu\] ■ V,,,n^ 

- 2 / PiV^,Aidtu)\ - dtu\) ■ V,,,(< - u\). 

Let us recall that the solitary wave r]2 satisfies the exponential decay estimate 

\d'^m{x -h- C2t)\ < c„e-'^(i+l^-'^-"2*l')^ for x,teR 
Consequently, we shall use the weight / defined by 

where e G [0, d] will be chosen sufficiently small. 
From ()7.2p . we first get the estimate 

\{[dt, Gm - Gi]u, u)\ < c(||/V,X/IIl2(5) + l|Vx,.« - 4)\\lHs) 

+ \\'^x,z{idtu)M - dtA^ - {dtu)\ + dtu\)\\L2^s^j ■ (\\Vx,zu\4\\L^i^S) + \Wx,zu\\\l^I^S) 

+ \Wx,z{{dtn)\,-dtu\,)\\L^s) + \\V^A{dtu)\-dtu\)U2^s))- (7.3) 



From the elliptic problem (j7.ip . we first get the estimates (this follows for example by using the 
decomposition (j3.6p of the solution, we refer to [J, [21] for example) 

l|V.,.n^||L2(5)<C|q3n|i2(ffi2)<C(|<pC/2|i2 + |C/2li2), i = M, 1, (7.4) 

where in the last step we used (|5.1ip . To conclude, we still need the estimates for fu\j, u\j — u\, 

{dtufj — dtu'j and {dtu)\,j — dtu\,^ — {dtu)\ + dtu\ with j = M, 1. We shaU deal with them one by 

one. 

1) Estimate for e-'(^+l^-''-''2ip)2 ^fe^^ ^^ -^ ^-^^ proof of Proposition^ we get that e-^(i+l^-''-'=2*P)2 ^^^ 

solves the elliptic equation 

V:,,, •PmV^.,,(/<^) = -[/, V^,^-Pa/V^.,J<^, in 5, 

/<U=o = fu, d^^'{fu\,)u=.i = -[/, arM]<,u=_i = 

We introduce the decomposition 

fu\, = m{z,\D\){fu)+v 

where m[z, \D\) is the Fourier multiplier ^[ '^1^ — —. We get for v the system 

'^x,z- Pm'^x,zV = -Vx,z- PM'^x,z{m{fu)) - [f,\/^^;,- PmVx,z]u\j, in S, 
v\z=o = 0, dzv\z=-i = 0. 
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From an energy estimate (as in the proof of Proposition 13. 2p . we obtain that 

+ (e'llM./ IIl2(5) + e\\V.AfuM)\\LHs))i\\v\\L^iS) + W"^ cc,zv\\lHS))) (7-5) 
where in particular we have use the fact that 

|V,,,/| < e/. 
Next, we can use the Poincare inequahty in the strip S which yields 

lblll,2{5) < C'I|Vx,2^'||l2(5) 
and the fact that u = Xif^2 which yields thanks to (jS.lip and Lemma [57 



\\miz, \D\){fu)h2^s) + l|V.,.(m(z, \D\)ifu))h2^s) < Ic{\^U2\l^ + \U2\l^] 
to obtain from (17. 5p by taking e small enough that 

\\V\\12 + \\V^,,V\\12 < ^C{\W2\l2 + \U2\l2). 



This yields 

2) Estimate for u\.j — u\. Let us set 



Vx,.(Mf)||i. + WfAAh < lc{\W2\l2 + \U2\I2). (7.6) 



g = -{Pm - Pi)V.,zU%, 

we get that u\^ — u\ solves 

V^,^ • PiV^,,«,^ - n5) = V^,, • g, in 5, 
u\i - u\ I ,=0 = 0, d, {uli -u\)\,=^i = Q 

The standard energy estimate for this problem yields 

||V,,,«, - u\)\\l2 < CMl < C{\\fuUl2 + ||v,.,,(/0||i.), 

which results (combined with the Poincare inequality and (|7.6p ) in 

IK - «?lli2 + iiv,,,(4, - U\)\\l2 < \C{\W2\12 + \U2\I2). 

3) Estimate for {dtu)^^ — dtvfj. As previously, we can set gj = dtPjVx,zu'j with j = M, 1, to get 
that [dtuf', — dtvf', satisfies the system 



yx,z-PjyxAidtu))-dtu)) = V:,^,-^j, in 5, 

{dtu)) - dtu%=o = 0, dl^{{dtu)] - dtu%=.i = -e, • g,|,=_i = 



and we obtain that 

\\{dtu)]-dtu)\\l2 + \\V,,.mu))-dtu])\\l2<CU,\^^^^ (7.7) 

thanks to fTI]) . 

4) Estimate for {dtu)\j - dtu\j - {dtu)\ + dtu\. We write v = {dtu)\,j - dtu\,^ - {dtu)\ + dtu\ 

here and we know from 3) that 

yx,z- Piyx,zV = Vx,z-[^M -^i- {Pm - PiWx,z{{dtu)\f- dtvll^)], in 5, 

^U=o = 0, dP^'vl=^i = -e^ • [gM - gi - {Pm - Pi)V x,z{{dtufM - dtu\j)]\z=-i 
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From a basic energy estimate, we obtain 

||v,,,^||2 < c(||(Pa/ - Pi)v.,,.{{dtu)l, - dtuii)\\L2 + \\dt{PM - Pi)v^,,uUy 

+ ||5iPiV,,,(<-n5)||i2 
< C[\\f{{dtu)i, - dtuiML^ + WV^MidtufM - dAML^ + \\f<\\L^ (7.8) 
+ ||V.,,(/OlU2 + ||V,,,(n^ - u'i)\\L2] 

where we use again / = e~'^(^+l^~'*~'^2i| )^ _ ^q conclude, we need to estimate w = f{{dtu)\j—dtu\^). 
We observe that w solves 

V,,, • Pa/V,,,u; = -[/, V,,, • PMV,Mdtu)ij - dtulj) - /V,,, • dtPM^x^zul^, in S, 
w\z=o = 0, dzw\z=-i = 0. 

We can proceed as in step 1) to obtain that 

\\w\\l, + \\Va:M\h < C\\fdtPMV.,Al\\h + \C{\W2\12 + \U2\I2) < ^C{\W2\l2 + \U2\I2) 



thanks to (j7.6|) . We can thus also obtain by combining the last estimate, (j7.8p and the estimates 
of step 1) and step 2) that 

Ibllia + \\V.M\h < ^^7(1^2112 + \U2\I2). 

Gathering all the previous estimates, we finally obtain (j5.12p . 

7.2. Proof of (I5J3D . Since 

{[[du Gm], Xi]u, v) = {[du Gm]xiu, v) - ([9j, Gm]u, xiv) 

= dt{GMXiu, v) - {GmXiu, dtv) - (Gudtixiu), v) - dt{xiGMU, v) 
+{{dtXi)GMU, v) + (xiGmu, dtv) + (xiGRidtu, v), 

we get by using the Green's Formula on the flat strip 5 that 

{[[dt, Gm], Xi]u, v) = dt [ PAiV^Miuf ■ yx,zv^ -dt f xi^A/V,,,n'' • V,,,vt 

Js Js 

-dt fiVx,zXi)-PM^x,zu''-v^- f PM^xAxiuf-VxAdtv)^ 
Js Js 

+ [ XiPaiV^^.u' • V^^Adtv)^ + / (V.,.xi) ■ PmV^.^^u'' {dtv)^ 
Js Js 

- [ PM^xAdtixiu))'' ■ V,,,t;t + / (5tXi)i^MV,,,n'' • V,,,t;t 
Js Js 

+ I Vx.z{dtXl) ■ Pm^x^zU^V^ + I XlPM^xAdtuf ■ Vx.zV^ 

Js Js 

+ I (yx,zXl)-PM'^xAdtufv^ 
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where u is u\,[ are satisfying ()7.ip with tjm, and u' = x(-2|-D|)7; with x a- smooth compactly 
supported cut-off function such that x(0) = 1- One can rewrite the above as 

([[dt, Gm], Xi]u, v) 

= I [dtPM^xAXluf ■ Vx,zV^ - dtPMXl^x,zu'' ■ Vx,zV^ - (V^,.Xi) • dtPMyx,zU^v^] 

+ / [PMy^Adtixiuf - (dtixiu))') ■ V.,zv^ - XiPM^xAidtu)' - dtu") ■ V,,,v^ ^^'^^ 



JS 

- (yx,zXi) ■ PM^xAidtuf - dtu') t;t] 
■.= Ai+A2 

while noticing that dtv' = (dtv)' . We will deal with Ai first. We have 

Al = [ dtPMVxMXluf - XlU^) ■ '^x,zV^ + / dtPMC^x,zXl)u^ ■ "^x.zV^ 

Js JS 

- j {yx,zXl)-dtPMVx,zU^v\ 

Since (xi^)^ ~ Xi^^ solves the equation 

yx,z ■ Pu^xAixiuf - Xiu^) = [Xi, V^,, • PmVxM^ in 5, 

{Xiuf - Xiu%=Q = dP'^'iixiuf - xiu')\z=-i = [Xi, dP'^']u\_^ = 0, 

by using the definition of xi and (j5.8p . we again obtain from an energy estimate that 

||V.,.((xin)'' - Xi^')IIl2 < \c{\^u\l2 + \u\l2). (7.10) 

By using again (15. Sp . we thus obtain 

Al < \c{\^u\l, + \^v\l, + \u\l, + \v\l,), 
where we used the fact that ||V2:^2i;T||2 < C|*Pi)|£^2. On the other hand, one can rewrite 
^2 = y PmVx,z [dtixiuf - {dt{xiu)f - Xidtu'' + Xi{dtuf) ■ V,,,t;t 

+ j Pm{Vx,zXi) ■ [{dtu' - {dtu)')Vx,zV^ - Vx,z{dtu' - (dtu)') v^ . 

We already have the estimate for dtu — {dtu) from step 3) of the last subsection, it only remains 
to estimate w = dt{xiuf' — {dt{xi'^))^ ~ Xidtu^ + Xi{dtu)^ ■ We get for w the equation in 5 

yx,z ■ PaNx,zW = - Vx,z • dtPMVxAxiu)^ + [Xi, ^x,z ' Pm"^ x,z]idtu^ - (dtu)'') 

+ Xl'^x,z ■ dtPM'^x,zU^ 

that we can rewrite as 

V^,^ • PmVx,zW = - Vx,z • dtPMyxAixiuf - Xiu'') + [xi, ^x,z • PMyx,z]{dtu'' - (dtu)^) 

+ [Xiyx,z-dtPMVx,z]u'' 

with the boundary conditions 

wU=o = 0, d^^'w\,=.i = 0. 
We thus get the estimate 

||t/;||i2 + \\Vx,zw\\l^ < c(\\VxAixiu)' - Xiu')\\l^ + ^Wdtu'' - {dtufWHi + IwAlm)- 
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Consequently, by using (j7.10p . ()7.7p and (j7.4p . we obtain 

\\w\\l2 + l|Vx,zU;||i2 < -(|<Pti|i2 + \u\l2) 
and hence 

|^2| < 1{Wu\l2 + |u|i2)(|*Pf|i2 + |7;|i2). 

Gathering the estimates for Ai and A2, we end the proof by substituting u = Xif^2 and v = U2 
into the estimates. D 
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